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ABSTRACT 

We present a systematic study of the physical properties of late-M dwarfs based 
on high-quality dynamical mass measurements and near-infrared (NIR) spectroscopy. 
We use astrometry from Keck natural and laser guide star adaptive optics imaging 



m 

(N) ; to determine orbits for the late-M binaries LP 349-25AB (M7.5-FM8), LHS 1901AB 

(M6.5-hM6.5), and Gl 569Bab (M8.5-FM9). We find that LP 349-25AB (Mtot = 

P^ I 0.120/qqq7 Mq) is a pair of young brown dwarfs for which Lyon and Tucson evolu- 

^^ ■ tionary models jointly predict an age of 140 it 30 Myr, consistent with the age of the 

>-Q ■ Pleiades. However, at least the primary component seems to defy the empirical Pleiades 

I . lithium depletion boundary, implying that the system is in fact older (if the parallax is 

^ . correct) and that evolutionary models underpredict the component luminosities for this 

^ , magnetically active binary. We find that LHS 1901AB is a pair of very low-mass stars 

(Mtot = 0-1941qq2^ -^0) with evolutionary model-derived ages consistent with the old 

age (> 6 Gyr) implied by its lack of activity. Our improved orbit for Gl 569Bab results 

^ ' in a higher mass for this binary (M tnt = O-l^O/gQgg Mq) compared to previous work 

0^ 



(0.125 lb 0.007 M0 lSimon et al.ll2006l ). We use these mass measurements along with our 
published results for 2MASS J2206-2047AB (M8-bM8) to test four sets of ultracool 
model atmospheres currently in use. Fitting these models to our NIR integrated-light 



o, 

f^ , dently from Dusty evolutionary models given the measured masses and luminosities. 



spectra provides temperature estimates warmer by ~250 K than those derived indepen- 
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We propose that model atmospheres are more likely to be the source of this discrep- 
ancy, as it would be difficult to explain a uniform temperature offset over such a wide 
range of masses, ages, and activ ity levels in the c o ntext of evolutionary models. This 
contrasts with the conclusion of iKonopackv et alj (|20ld ) that model-predicted masses 
(given input T^s and Lboi) are at fault for differences between theory and observations. 
In addition, we find an opposite (and smaller) mass discrepancy from what they report 
when we adopt their model-testing approach: masses are too high rather than too low 
because our T^s estimates derived from fitting NIR spectra are ~650 K higher than 
their values from fitting broadband photometry alone. 



Subject headings: binaries: close — binaries: general — binaries: visual — infrared: 
stars — stars: low-mass, brown dwarfs — techniques: high angular resolution 



Introduction 



While the spectra of late-M dwarfs are by definition very similar, their underlying nature can 
vary widely from youthful brown dwarf s to old low- mas s stars . The first field dwarf with a spectral 



type later than M6 was discovered by Ivan Biesbroeckl (11944 ), yet even 50 years later the nature 



of such objects was open to debate, with iKirkpatrick et al. 



199J) suggesting that the latest type 



M dwarfs may all be young enough to be substellar. iRebolo et al.l (jl992l ) proposed a method for 
discriminating between brown dwarfs and low-mass stars by using the Li I doublet at 6708 A, 
as models predict that stars and the highest mass brown dwarfs (M > 0.06 Mq) should deplete 
their initial lithium rapidly. Such model predictions of lithium depletion have al so been used to 
age-date open clusters (e.g.. IStauffer et al.lll998l : iBarrado y Navascues et al.lll999l ). However, the 
underlying theory still remains unconstrained by direct mass measurements for ultracool dwarfs 
with lithium measurements near the predicted depletion boundary. In recent years, spectroscopic 
signatures of low surface gravity have been used to ident ify very young (-^y 10-100 Myr) ultracool 
dwarfs among the field population of late -M and L dwarfs (JMcGovern et al.ll2004l : lAllers et al.ll2007l : 
Kirkpatrick et al.ll2008l : iRice et al.ll2010l ). thus enabling a means of discriminating between stars 
and brown dwarfs if the objects are young enough. 

Dynamical mass measurements provide a direct means for studying the diversity in the phys- 
ical properties of late-M dwarfs. We present here masses for the late-M binaries LP 349-25AB, 
LHS 1901AB, and Gl 569Bab based on Keck natural guide star (NGS) and laser guide star (LGS) 
adaptive optics (AG) imaging from our ongoing orbital monitoring program targeting ultracool 
binaries. We combine our new measurements with published results to perform a systematic study 
of late-M dwarfs. To that end, we also present medium-resolution (R ~ 2000) near-infrared (NIR) 
spectra for this sample of late-M binaries. The resulting combination of high-quality dynamical 
masses and spectra enables strong tests of model atmospheres over a wide range of masses, ages, 
and activity levels. 



LP 349-25 was first identified as a nearby M8 dwarf bv lGizis et al.1 (l200d). andlForveille et al 



( 2005 ) discovered its binary nature after being initially unresolved bv lClose et al.l (J2003l ). iGatewood k. Coban 



2009 ) measured a trigonometric parallax of 75.8 it 1.6 mas for LP 349-25, corresponding to a 



distance of 13 .2 ifc 0.3 pc, which is somewhat more distant than the photometric estimate of 
10.1 ± 1.2 pc (jForveille et al.l l2005l ) . iReiners &: Basril (J2009l ) obtained high-resolution integrated- 
light spectroscopy using Keck/HIRES and found no evidence for lithium (equivalent width < 0.5 A). 
Konopacky et al.l (120101 ) recently published a dynamical mass of 0.122 it 0.009 Mq based on Keck 



LGS AO data. We derive here a dynamical mass of 0.120;to'oo7 ^©) based on an independent set 
of astrometry. Thus, this binary is a pair of brown dwarfs. 



LHS 1901 was first identified as a nearby M 6.5 dwarf bvlReid et al 



([20091) have independently classified it as M7.0. iMontagnier et al. 



(|2003l ). and iLepine et al 



(120061) discovered LH S 1901 to 



be a 0'.'2 binary, which was also later identified by lLaw et al.l ()2008l ). iLepine et al.l ()2009l ) measured 
the parallax of LHS 1901 to be 77.8 ± 3.0 mas, corresponding to a distance of 12.9 ± 0.5 pc. Our 
dynamical mass of 0.194^q'q2;^ Mq indicates that this is a pair of very low-mass stars in a highly 
eccentric orbit (e = 0.830 ± 0.005). 

Gl 569B was disc overed as a companion to the chromospherically active M2 star Gl 569 
by iForrest et al.l ( 19881 ). and Martin et al.l (J2000l ) resolved it as an M8.5+M9 binary using Keck 



AO. There are three published orbit det erminations for this binary (IZapatero Osorio et al. 



2004 



Simon et al.ll2006l : iKonopacky et al.ll2010l ) that all give consistent values of the total mass (0.125 it 
0.007 Mq). However, these orbits are all based largely on astrometry from 1999 to 2001 obtained 
with the Keck NIR cameras KG AM and SCAM (the slit-viewing camera for NIRSPEC). Since these 
instruments are not astrometrically well-calibrated, we have derived a new orbit for Gl 569Bab 
based solely on data from Keck/NIRC2 and the Hubble Space Telescope (HST) Space Telescope 
Imaging Spectrograph (STIS). In addition , we use the revised Hipparcos parallax for Gl 569A 
of 103.59 ± 1.72 ma s (Ivan LeeuwenI 120071 ) rather than the original value of 101.91 it 1.67 mas 



(jPerrvman fc ESAl 119971 ) used in all previous orbit determinations. Our improved measurement 



of the total dynamical mass for Gl 569Bab is 0.140_QgQo Mq, indicating that the binary is more 



massive than previously thought but still with at least one unambiguously substellar component. 

In addition to these three binaries, we include the M8+M8 binary 2MASS J22062280-2047058AB 
(hereinafter 2MASS J2206-2047AB) in our study as it also has a well-determined mass (0.15tQ;Q| M, 
Dupuv et al.ll2009al ) and integrated- light spectroscopyo 



0) 



^The only other late-M dwarf binary with a precise dynamical mass is LHS 1070BC (A/tot ~ 0.157 ± 0.009 Mq 
Leinert et aUboOll : [Seifahrt et allboosl ). but it is only I'.'S from the bright primary LHS 1070A (M5.5) making it 
difficult to obtain spect roscopy of the late-M dwarfs. We also do not consider four other late-M dwarf binaries from 
Konopacky et al.1 (2010|), as their dynamical masses are poorly constrained with 60-190% errors. 



Observations 



2.1. Keck/NIRC2 AO 

We have monitored LP 349-25AB, LHS 1901AB, and Gl 569Bab with the AO system at the 
Keck II Telescope on Mauna Kea, Hawaii, using the facihty NIR camera NIRC2 in its narrow field- 
of-view mode. At each epoch, we obtained data in one or more filters coverin g the standard atmo- 
spheric windows from the M auna Kea Observatories (MKO) filter consortium (ISimqns &: Tokunaga 

2002 : lTokunaga et al.l[20oi '). For LP 349-25 AB, we used LGS AO (|wizinowich et al.ll20od : lvan Dam eTal 
20061 ) when observing conditions allowed, and NGS AO otherwise (see Tabled]). The LGS bright- 
ness, as measured by the flux incident on the AO wavefront sensor, was equivalent to a V ~ 9.5- 
10.4 mag star. The tip-tilt correction and quasi-static changes in the image of the LGS as seen 
by the wavefront sensor were measured contemporaneously by a second, lower bandwidth wave- 
front sensor monitoring LP 349-2 5, which saw the equiva l ent ol an R rj 14.3-1 4.6 mag star. For 



LHS 1901 AB, we used NGS AO (IWizinowich et al. 



200C 



van Dam et al.ll2004l ). and the incident 



flux on the wavefront sensor from LHS 1901 was equivalent to a y ~ 13.0-13.4 mag star. For 
Gl 569Bab, the primary star Gl 569A (R = 9.4 mag) located S'.'O from the science target provided 
the NGS AO correction. 

Our procedure for obtaining , reduc ing, and analyzing our images is described in detail in our 



previous work (e.g.. lDupuv et al.ll2009bH d). Table [T] summarizes our observations of LP 349-25AB 



LHS 1901 AB, and Gl 569Bab. The binary separation, position angle (PA), and flux ratio were 
determined using the same three-component Ga ussian representation of the point-spread func- 
tion (PSF) as described in iDupuv et al.l (J2009d ). At epochs where the binary was sufficien tly 
well-separated, we were also able to use the StarFinder software package (IDiolaiti et al.ll2000l) t o 
simultaneously solve for the PSF and binary parameters. As described in iDupuy et al.l (|2009d ). 
we assessed systematic errors in both PSF-fitting procedures by applying them to artificial bi- 
nary images constructed from images of PSF reference stars with similar FWHM and Strehl. 
We found good agreement between both methods, with S tarFinder typicall y giving equivalent 
or smaller errors. We adopted the astrometric calibration of lGhez et al.l (J2008l ). with a pixel scale 
of 9.963ib0.005 mas pixel"^ and an orientation for the detector's +y-axis of +0?13ib0?02 east of 
north. We applied the distortion correction developed by B. Cameron (2007, private communica- 
tion), which changed our astrometry below the la level. The resulting relative astrometry and flux 
ratios for LP 349-25AB, LHS 1901AB, and Gl 569Bab are given in Tabled 

For Gl 569Bab, we also used our Keck images to measure relative photometry between the 
primary star Gl 569A and Gl 569Bab at JHK, as there is no photometry published in the MKO 
system. We summed the flux in circular apertures with radii of 0'.'3 to 0'.'5. We used NIRC2 in 
subarray mode to reduce the minimum allowed exposure time and thus avoid saturating the primary, 
which limited the size of the largest aperture we could use. However, in experimenting with even 
smaller apertures we found the same relative photometry to within la. We used a portion of the 
array minimally affected by light from either the bright primary star or the binary to measure the 



median sky level, which was subtracted from our summed aperture fluxes. The relative flux of the 
binary in integrated light compared to Gl 569A was A J = 4.15 it 0.05 mag, AH = 4. 14 it 0.03 mag, 
and AK = 3.86 it 0.03 mag. To compute absolute photometry for Gl 569B we applied these flux 
ratios to MKO photometry for Gl 569 A, which was derived by converting its Two Micron All Sky 
Survey (2MASS) photometry to the MKO system using offsets computed from synthetic photometry 
of a SpeX prism spectrum of Gl 569A (A. Burgasser 2010, private communication. The conversions 
from 2MASS to MKO for J, H, and K were —0.03 mag, +0.02 mag, and —0.01 mag, respectively 
(all < 1.3(7 compared to the 2MASS errors). Thus, the final MKO integrated-light photometry for 
Gl 569B was J = 10.75 ± 0.06 mag, H = 10.15 ± 0.04 mag, K = 9.62 ± 0.03 mag. 



2.2. HST/STIS Imaging 

Gl 569Bab was observed by HST/STIS on 2002 June 26 UT as part of a program to obtain 
resolved optical spectroscopy of ultracool binaries (GO-9499; PI Martin). During acquisition, the 
STIS CCD imager took two images with the long-pass filter F28X50LP in which the binary is 
well-detected and the primary star Gl 569A is saturatedo We anal yzed these images using our 
PSF-fitting routine base d on the TinyT i m model of the HST PSF (lKristill995l ). as described in 
our previous work (e.g., iLiu et al.l l2008l : iDupuy et al.l l2009al ). We adopted the pixel scale and 
distortion solution from the STIS Instrument Science Report 2001-02 (jwalsh et ahlboOll . 50.725 ± 
0.056 mas pixel"^). 



We tested for systematic errors in our best-fit binary parameters by applying our PSF-fitting 
routine to artificial binary images created from STIS CCD images of single stars from HST cali- 
bration programs in Cycles 8 and 10 (CAL/STIS-8422 and 8924). To preserve the undersampled 
nature of the PSF, we created only artificial binaries separated by an integer number of pixels. 
The separation and instrumental PA of Gl 569Bab in the STIS images were 1.92 pixels and 3?1, so 
artificial binaries with A{x,y) = (0,2) were an excellent representation of the science data. The 
single stars were ~10x brighter in the STIS images than Gl 569Ba, so we scaled them to match the 
science data and then injected them into the actual science image to accurately simulate the noise. 
The binary images were injected at a distance of 5" from Gl 569A, chosen to be comparable to 
Gl 569B's separation. For the final adopted errors we used the larger of the rms from the artificial 
binaries or the rms of the best-fit parameters to the science data. From our artificial binary tests, 
we found systematic offsets in separation (3.0 mas, 2.3cr), PA (0?2, 0.2cr) and fiux ratio (0.15 mag, 
5a), which we applied to our best-fit binary parameters for Gl 569Bab. The final separation, PA, 
and flux ratio we derived were 97.3 ± 1.3 mas, 94?0 ± 0?9, and 0.99 ± 0.03 mag (Table E]). 



F28X50LP has a cut-on wavelength around 0.55 /im and extends as red as 1.0 fj,m. 



2.3. VLT/NACO 

We retrieved archival images of LP 349-25AB obtained with the Very Large Telescope (VLT) 
at Paranal Observatory on five epoc hs spanning six mon t hs in 2006. These d ata were taken with 



the NACO adaptive optics system (jLenzen et al. 



2003 



Rousset et alJ l2003l l using the N90C10 



dichroic and S13 camera. We adopt the pixel scale of 13.221 it 0.017 mas pixel"^ derived by B. 
Sicardy (2008, private communicationjfl from observations o f Pluto's motion coni pared to its known 
ephemeris, and an orientation of the detector of 0?0 ± 0?1 (JForveille et al.ll2005l ). We followed the 



same procedure for analyzing these data as in our previous work with VLT images (JDupuy et al 
2009al Jd). We registered, sky-subtracted, and performed cosmic ray rejection on the raw archival 



images. We used the same analytic PSF-fitting routine as was used for the Keck data to fit the 
VLT images, with uncertainties determined from the standard deviation of measurements from 
individual dithers. The resulting best-fit parameters are summarized in Table [2l 



2.4. Previously Published Astrometry 



Forveille et al.l (120051 ) published two epochs of astrometry for LP 349-25AB, and we investi- 
gated the possibility of improving the precision of their measurements by applying our PSF-fitting 
procedure to their data, which is in public archives. The CFHT observations of LP 349-25AB on 
2004 July 3 UT were not available in the CFHT archive because they were obtained during engi- 
neering time (T. Forveille 2010, private communication), but the VLT/NACO observations from 
2004 September 24 UT are in the VLT archive. We measured a separation of 111.2 it 0.3 mas, a PA 
of 5.66 ± 0.12°, and a flux ratio at H ba nd of 0.40 =fc 0.03 mag. Our flux ratio and separation are in 
good agreement with those measured bv lForveille et al.l (|2005l ). 0.38 it 0.05 mag and 107 it 10 mas, 
while our PA of 5.66 it 0.12° is somewhat inconsistent with their 7.1 ± 0.5° (at 2.8o"). We tried 
using each of these measurements to fit the orbit and found that they gave consistent dynamical 
masses with a very similar x^ for the best-fit. Our higher precision measurement would enable the 
uncertainty in some of the orbital parameters to be imp roved by up to a facto r of ~2. However, we 
conservatively chose to use the published astrometry of lForveille et al.l (120051 ) as their larger errors 
and PA offset may result from analysis of astrometric calibration data unavailable to us that was 
taken along with their observations. 



Montagnier et al.l (J2006l ) published three epochs of astrometry for LHS 1901 AB, and the first 



two data sets are available in the CFHT archive. Unfortunately, not all of the astrometric calibration 
data was available in the archive, so independent measurements from these data cannot be used in 
the orbit fit becau se, e.g., the orientation of the detector at the time of those observations is not 
accurately known. iMontagnier et al.l (J2006l ) found uncertainties of 5 mas and 0.5° for all three data 
sets. The rms scatter of our PSF-fitting measurements is comparable or slightly lower (3-4 mas 



Ihttp : //www . eso . org/sci/f acilities/paranal/instruments/naco/doc/VLT-MAN-ESD- 14200-2761_v83 . 3 . pdf | 



and 0.4-0.6°), and thus even with the needed calibratio n data a drastic improve ment in precision 
seems unhkely. Thus, we use the pubhshed astrometry of lMontagnier et al.l (|2006l ) when fitting the 
orbit of LHS 1901 AB (Table [2]). 

For Gl 569Bab, we do not use any of the previousl y published a s trome try. We performed our 
own analysis on the ori ginal Keck/NIRC2 imag es from ISimon et al.l ( 20061 ) following the method 
described in Section [2.11 iKonopackv et al.l (J20ld ) published one epoch of astrometry for Gl 569Bab 
which is contemporaneous with our 2009 May data; we conservatively choose to exclude their data 
from our orbit fit, although it is consistent with ou r best-fit o r bit. W e do not include any of 



the K eck/KCAM or SCAM astrometry published by iLane et al.l (|200ll ) or IZapatero Osorio et al 



(J2004l ). as those cameras have not been rigorously astrometrically calibrated, unlike NIRC2 and 
the STIS CCD. 



2.5. IRTF/SpeX Spectroscopy 

In addition to astrometry, which we use to determine orbits and thus dynamical masses, we have 
also obtained NIR spectra to perform tests of model atmospheres, as discu ssed in Sections [3] a. nd [H 
Using the NASA Infrared Telescope Facility (IRTF) spectrograph SpeX (JRavner et al.l 120031 ) . we 
obtained spectra of LP 349-25, LHS 1901, and Gl 569B on 2008 July 6, 2010 January 30, and 2010 
March 3 UT, respectively. We used SpeX in cross-dispersed (SXD) mode, with five orders spanning 
0.81-2.42 /im. For LP 349-25 we used a 0'.'5 slit (R = 1200), and for LHS 1901 and Gl 569B we used 
a 0'.'3 slit {R = 200 0). We calibrated , extracted, and tellur ic-corrected the data using the SpeXtool 
software package (IVacca et al.l 120031 : ICushing et al.l 1200J) . We used integrated- light photometry 
from 2MASS to flux-calibrate our spectra of LP 349-25 and LHS 1901. For Gl 569B, we used our 
own integrated-light photometry derived from our Keck images (Section 12. ip to flux calibrate its 
spectrum. We used the mean of the JHK scaling factors computed for each bandpass to place our 
spectra on an absolute flux scale. 



3. Results 



3.1. Orbit Determination and Dynamical Masses 

Combining our Keck AG monitoring with published discovery data, we have observations span- 
ning 5.9 years for LP 349-25AB, 5.9 years for LHS 1901AB, and 7.9 years for Gl 569Bab. Following 
our previous studies, we have determined the orbital parameters and thei r uncertainties us ing both 
a Markov Chain Monte Carlo (MCMC) approac h, described in detai l by lLiu et al.l (120081 ) . as well 
as the least-squares minimization routine ORBIT (JForveille et al.lll999l ). In addition, we have devel- 
oped our own visual binary orbit fitting routine that utilizes the Levenberg-Marquardt algorithm 
as implemented in the MPFIT routine for IDL (JMarkwardtl l2009l ) . To determine uncertainties on 



the best-fit parameters from this routine, we have performed Monte Carlo simulations in which 
Gaussian noise is added to the data corresponding to their measurement uncertainties. We used 
the resulting distributions to determine the medians and standard deviations of the best-fit pa- 
rameters. When multiple astrometric measurements were available from different bandpasses at a 
single epoch, we used the highest precision measurement in the orbit fit (see Table [2]). The best-fit 
orbits are shown in Figures [T] and EJ 

We found essentially symmetric uncertainties in the MCMC-derived orbital parameters for 
LP 349-25AB. These parameters agreed to within < 0.5% of those derived by our custom fitting 
routine, and the errors agreed to within < 15%. As an additional test of our Levenberg-Marquardt 
fitting routine, we found that the best-fit parameters from ORBIT also agreed to within < 1%, and 
the errors derived from our Monte Carlo approach agreed to within < 10% of the errors computed 
from the full covariance matrix in ORBIT. The orbital parameters and errors derived by our MCMC 
analysis and custom fitting routine are given in Table [3l 



We experimented with including the five epochs of relative astrometry from iKonopackv et al 



( 20ld ) when fitting the orbit of LP 349-25 AB. If all of their data are included, the x^ becomes 
unrealistically large: x^ of 58.2 for 31 degrees of freedom (DOF). This large x^ value is due to 
the fact that two of their measurements from 2007 December and 2008 May are extreme out- 
liers from the best-fit orbit. (Indeed, their own orbit fit is inconsistent with these data as the 
Konopacky et al.ll2010l orbit has a reduced x^ of 2.15, which is unrealistically large.) By excluding 



their most discrepant epoch (2008 May 30 UT), we were able to achieve a reasonable x^ (33.1 for 
29 pOF). Neither the orbital parameters nor their uncertainties change significantly after including 



the IKonopackv et al.l (|20ld ) data. This is because their astrometry does not offer additional unique 
constraints on the orbit, as it was obtained contemporaneously with o urs and is of lowe r preci sion 



by a factor of ~4 to 10. Thus, we conservatively choose to exclude the IKonopackv et al.l (J2010l ) as- 
trometry from our analysis. Our best-fit total mass (Mtot) for LP 349-25AB is 0.1205 ±0.0007 M© 
(0.6%), and after including the 2.1% error in the parallax this becomes 0.120J:o:oo7 M© (6%). 

From our MCMC analysis of LHS 1901 AB, we found significant asymmetries in the posterior 
probability density distributions of the orbital period, semimajor axis, and eccentricity. The con- 
fidence limits on all orbital parameters are given in Table [3l and the distributions of these three 
parameters are shown in Figure [3l The orbital pe riod and sernimaior axis display a high degree 



of covariance as expected (e.g., see Section 3.2 of iDupuv et al.l l2009bl ) . There appear to be two 
slightly degenerate solutions, with one having a shorter period, lower eccentricity, and higher mass 
than the other. Our MCMC analysis enables us to treat this degeneracy properly by accounting 
for it in the confidence limits of these orbital parameters and the resulting dynamical mass. The 
difference in mass between these two solutions is very small (about 1.6%), insignificant compared 
to the 12% error contributed by the distance uncertainty. Our best-fit total mass for LHS 1901AB 
is 0.1944 lb 0.0028 Mq (1.4%), and after including the IL^'iVo error in the parallax this becomes 
0.194tO:J]25 Mq (12%). 



The orbital parameters we find from Gl 569Bab are all consistent with previously published 
orbits, with a singular exception. T he semimajo r axis w e find (95 .6l^ 'n mas) is 4.3c7 higher tha n the 
value of 90.4 ± 0. 7 mas derived bv ISimon et al.1 (J2006l ). and the IZapatero Osorio et al.l (J2004l ) and 
Konopackv et alj (J2010|) values are similarly discrepant (by 1.6a and 3.7a, respectively). Since 
all orbit determinations to date have been dominated by the nine epochs of astrometry from 
Keck KCAM and SCAM, the semimajor axis discrepancy implies that the pixel scale of one or 
both of these instruments was underestimated by ~5%, resulting in systematically low semima- 
jor axes. This is supp orted by the fact that the previously published NIRC2 astrometry from 



Simon et al.l (|2006l ) and lKonopackv et al.l ((20101) is in good agreement with our best-fit orbit. Also, 
the HST/STIS measurement would be highly inconsistent with a value for the semimajor axis as 
low as previously found, showing that our semimajor axis is not simply due to a problem with the 
NIRC2 calibrationo Since we find the same orbital period as in previous work but a semimajor 
axis that is 5.4% larger, our dynamical total mass is significantly higher than previously found. 
This would be a 16% effect, but it is moderated somewhat by the revised Hipparcos parallax that 
is 1.6% (about la) higher than the original value, which brings the mass down by 4.8% for a net 
change of 11%. Our best-fit total mass for Gl 569Bab is 0.140]'^q'qq| Mq, and after accounting for 
the 1.6% parallax error this becomes 0.140^'^Q'gQg Mq. Our MCMC analysis also reveals a slightly 
degenerate orbit solution with double-peaked distributions for the argument of periastron (a;) and 
PA of the ascending node (Jl) as shown in Figure [3l The degeneracy in these values are apparently 
not correlated with any other orbital parameters. 



3.2. Spectral Types 



Gizis et al.l (|200d ) used the integrated-light optical spectrum of LP 349-25 to derive a spectral 



type (SpT) of M8.0 it 0.5. Comparing our near- infrared spectrum o f LP 349-25 to spec tra of M 



dwarf standards from the IRTF Spectral Library (jCushing et al. 



2005 



Rayner et al.ll2009l ). we find 



an infrared spectral type of MB. Without resolved spectroscopy of the binary, we cannot directly 
determine the spectral types of the individual components. However, we have used our photom- 
etry and the emp i rical S pT-absolute magnitude relations of I Cruz et al.l (120031 . J band only) and 



Golimowski et al. 



(120041 . L' band only) to estimate spectral types. We used the method described 



bv iDupuv et al.l (|2009d ). which accounts for both the measurement uncertainties and the intrinsic 
scatter in the empirical absolute magnitude relations. We found that the ra nge of J-band abso lute 
magnitudes allowed by our measurement errors exceeded the bounds of the ICruz et al.l (J2003l ) re- 
lation, and thus we base our estimates only on the L'-band relation. We estimate spectral types of 
M7.5ib 1.0 for the primary and MS.Olio ^°^ ^^^ secondary, which are consistent with the published 
integrated-light spectral type. 



*We also note that a calibration problem with NIRC2 would have been apparent in our previously published binary 
orbits, which typically have data from NIRC2 as well as other well-calibrated instruments such as ifST/WFPC2-PC, 
HST/ ACS-HRC, and VLT/NACO. 
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Reid et al.1 (120031 ) measured the integrated-light spectral type of LHS 1901 to be M6.5, and 
Lepine et alj ( 2003 ) measured M7.0. Our near-infrared spectrum is better matched by the M dwarf 
standard vB 8 (M7) than GJ 1111 (M6.5), thus we find an infrared type of M7 for LHS 1901. 
Using the method described above, the J-band absolute magnitudes of the individual components 
provide an estimated spectral type of MG.Slg'g for both the primary and secondary, consistent with 
measured integrated-light spectral types. The difference of zero subtype between the components 
is also consistent with our measured flux ratios and late-M dwarfs with known parallaxes. For 
example, the rms scatter among i^-band absolute magnitudes for M6.5 dwarfs is 0.12 mag (cf. 
AK = 0.107 ± 0.007 mag for LHS 1901AB). 



For Gl 569Bab, we adopt the spectral types of M8.5±0.5 and M9.0±0.5 derived bv lLane et al 



(|200ll ). They used their resolved near- infrared spectra of the individual components to compute 
spectral indices and compared these values to field dwarfs of known spectral type. They found a 
difference of 0.5 s ubtypes between the two com ponents and assigned the integrated-light spectral 
type of MB. 5 from lHenry &: Kirkpatrickl (jl990l ) to the primary, resu lting in a type of M 9.0 for the 
secondary. These spectral types are consistent with those found by lMartin et al.l (120061 . M9+M9) 
using resolved optical spectroscopy of Gl 569Bab from HST/STIS. 



3.3. Bolometric Luminosities 



We computed the total bolometric luminosities (Lboi) of our sample binaries from their integrated- 
light spec tra and photometry . In ad dition to our SpeX spectra, we used the L' and M' band p h otorn - 
etry from iGolimowski et al.l (|2004l ) and the 24 //m Spitzer photometry fromlGautier et al.l (120071 ) 
for LP 349-25, and we used the L' band photometry from iForrest et al.l (jl988l ) for Gl 569B. For 



LHS 19 01, we estimated a K — L' color of 0.56ib0.07 mag from the empirical relation of lDupuy et al 



(|2009al ) to extend its SED to longer wavelengths. To accurately account for the SED at shorter 
wavelengths, we appended optical spectra, using the overlapping region of our SpeX spectra to 



determine the flux scaling. For LP 349 -25 and Gl 569B we used optical spectra from iReid et al 



(|2003l ) and iHenrv &: Kirkpatrickl (Il99d). respectively. For LHS 1901 we used the optical spectrum 
of the M7 standard vB 8 (JHenrv &: Kirkpatricklll990l ). We numerically integrated the spectra and 
photometry points, interpolating between the gaps in the data, extrapolating to zero flux at zero 
wavelength, and assumi ng a Rayleigh-Jeans t ail beyond the last photometric point. We accounted 
for errors as described in lDupuy et al.l (l2009al ). finding a total log(Lboi/-^0) of - 
LP 349-25, -2.67+^-^t dex for LHS 1901, and -3.212 ± 0.018 dex for Gl 569B. 



2.804+°:°^6 dex for 



We apportioned these total bolometric fluxes to the individual components by using the mea- 



sured K-hand flux ratio and the bolometric correction-SpT relation derived by IGolimowski et al 



(| 20041 ). properly accounting for uncertainties in the component spectral type estimates. The re- 
sulting individual luminosities are given in Table HI For the componen ts of Gl 569Bab, our L^oi 
values are in good agreement with the estimates from lLane et al.l (J200ll ). As in our previous work, 
we correctly account for the covariance between the individual luminosities (via the flux ratio) in 
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the following analysis, allowing more precise determinations of relative quantities (e.g., AT^fj and 
the mass ratio). 



3.4. Tgfj and log{g) from Model Atmosphere Fitting 

We used a variety of solar-metallicity model atmosphere grids to fit our integrated-light 
IRTF/ SpeX spectra, , sum marized in Table [5j Th e two publicly available g rids we used are Ames- 
Dusty (JAllard et ahlboOlh and PHOENIX-Gaia (JBrott k Hauschildtl I2OO5I ) . We also used an im- 
proved version of Ames-Dusty that employs the same treatment of dust but with u pdated line 
lists for some important opacity sources such as FeH, CrH, and H2O (see Section 3.2 o f 



Rice et al. 



20ld ): this "Gaia-Dusty" model grid is the same as that employed bv lKonopackv et al.l ([20101). Our 
fourth model grid Drift-PHOENIX (jWitte et al.ll2009l ) relies on yet another gas opacity data set 
and, more importantly, features improved dust opacities. The model treats the dynamics of the 
dust in a presently unrivaled level of detail by calculating the rate of seed format ion and solving an 
actual growth rate equation system for gravitatio nal settling of composite grains (jWoitke &: Helling 
20031 : iHelling fc Woitkj [200^ : iHelhng et allboosi ). 



We fit model atmospheres to our spectra using the procedure outlined by iBowler et al.l (J2009l ) . 
For every model spectrum, we found the optimal scaling factor to match the observed spectrum 
by minimizing the x^ statistic. The model spectrum with the lowest resulting x^ value gives the 
best-fit parameters for Tgfj and log{g). For each spectrum, we fit individual wavelength ranges 
corresponding to the standard near- infrared bandpasses (YJHK). We also performed fits to the 
entire SED, both with and without the segment of the spectrum between 0.81 and 0.95 um. This 



regio n of the spectrum is often excluded when fitting models to L and T dwarfs (e.g.. lCushing et al 



20081 ) because it is strongly affected by the broad wings of the resonant K I doublet at 0.77 /im, 
which have lon g been difficult to rn odel (e.g.. lBurrows &: Volobuyevll2003l ) and can be impacted by 
dust modeling (jjohnas et al.ll2008l ). Such alkali lines are present but weaker in late-M dwarfs, and 
they may cont ribute to observed d iscrepancies with models between ~0.7 and 0.9 /im for objects 
as early as LO (JReiners et al.l 120071 ). 



The results of our model atmosphere fitting procedure are given in Table [H We report fits 
both for our full wavelength range (0.81-2.42 //m) and for the NIR only (0.95-2.42 //m). We prefer 
the NIR SED fits for the following reasons: (1) the NIR SED fits are more consistent with the 
individual YJHK band fits (typically within one model grid step of 100 K, 0.5 dex in log{g)); 
(2) the full SED fits are systematically cooler by ~100 K compared to the individual band fits, 
implying that the lowest S/N portion of the spectrum from 0.81 to 0.95 /im has a disproportionate 
influence on the results; and (3) by adopting a 0.95 /im cutoff here, our results will be directly 
comparable to future work on L and T dwarfs that require such a cutoff because of the problems 
with K 1 lines. 



Figure H] shows the x^ contours of the flts to our observed NIR spectral energy distributions 
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(SEDs) using all model spectra. The elongation of the contours in the log(g') direction indicates that 
gravity is less well-constrained by our fits than effective temperature. We also note that there are 
occasionally local ^ minima in these plots that correspond to physically implausible parameters. 
For Drift-PHOENIX, the deviant x^ minimum occurs at very low temperatures (2000 K) and high 
gravities (log (5) = 6.0, which for masses of <0.08 Mq correspond to radii of <0.045 i?Q, more than 
a factor of two lower than evolutionary model predicted radii). For the other three model grids, 
~)(^ minima extend to higher temperatures (~3300 K) and extremely low gravities (log((7) = 3.5, 
which for masses of >0.06 Mq correspond to radii of >0.7 i?©). In the rare cases when the global 
minimum is found in one of these physically implausible regions (~3% of fits, see Table ED, we 
select the second best fitting model instead, and this always results in best-fit parameters from the 
prevailing solutions of Tefj = 2700-3000 K. 

The different grids of model atmospheres generally give very similar best-fit parameters despite 
their significantly different input assumptions. For example, the treatment of dust in the Drift- 
PHOENIX models is much more realistic than in other models, yet for our sample it yields the 
same Tefj values as other grids. Thus, it appears that the effects of dust at the warm temperatures 
of late-M dwarfs only subtly impact the resultant model spectra for our purposes. In fact, for the 
NIR SED fits (0.95-2.42 //m), all four sets of models give exactly the same Teg. 

Finally, we investigated how binarity impacts the atmospheric parameters derived from integrated- 
light spectra. We constructed artificial binary integrated- light spectra by summing individual 
IRTF/SpeX SXD spectra, both from our own spectra (e.g., summing the spectra of LP 349-25 and 
Gl 569B) and from p ublicly available spe c tra of the M7, M8 , and M9 spectral standards vB 8, 



vB 10, and LHS 2924 (ICushing et al.ll2005l : iRavner et al.ll2009l ). We fit these spectra both individ- 



ually and summed, after scaling down the later type spectrum to best match the measured JHK 
flux ratios of LP 349-25AB and Gl 569Bab|j Comparing the summed and input "primary" spec- 
trum fits, we found that the best-fit values from individual bands were nearly identical (at most 
different by one grid step of 100 K or 0.5 dex). In some cases of we found that the summed spectra 
gave systematically lower temperatures of 100 K as compared to the input "primary" spectrum. 
However, this was not the case when the input spectra had JHK colors comparable to the compo- 
nents of LP 349-25 AB and Gl 569Bab. We note that within each binary the components have very 
similar JHK colors, so we might expect that simulated component and summed spectra would give 
the same results. For example, the integrated-light J — K color of LP 349-25 is 1.02 mag, and its 
components are only different by ~0.03 mag {J — K = 1.00 mag and 1.05 mag) despite the binary 
fiux ratio of ~0.3 mag (the H — K and J — H colors are even closer, Acolor ~ 0.01 mag). 

In summary, our simulations show that the individual band fits should be unaffected by con- 



^LHS 1901 and 2MASS J2206-2047 both have flux ratios within f«0.1 mag of unity across JHKL' bandpasses 
and thus have components that are likely to be spectrally identical. For example, LHS 1901 has a bolometric flux 
ratio of 0.044 ± 0.015 dex, and thus the difference in T^b between the individual components is expected to be 75 K 
(i.e., less than one model grid step) for an assumed radius ratio of unity. 
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tamination from the spectrum of the secondaries in our sample binaries, and even if the SED fits 
are affected they will give Tgfj values systematically lower by only 100 K. Given that our best-fit 
values of the NIR SEDs show no such systematic offsets from the band fits, we conclude that our fits 
to the integrated-light spectra accurately find the best-fit parameters of the primary component's 
spectrum. 



3.5. Age Constraints from Kinematics and Activity 

We have considered whether the space motion and activity of our sample binaries can pro- 
vide useful constraints on their ages. For L P 349-25, there a r e two discrepant measurements of 
its radial velocity (RV) in the literature: (1) iReiners Sz Basril ( 20091 ) derived —16.8 ± 2.0 km s~^ 
fro m cross-correlation of t heir Keck/HIRES spectrum with the RV standard Gl 406 (M6); and 
(2) iKonopackv et al.l (2010|) found a center-of-mass velocity of —8.0 ±0.5 km s~^ by fitting an orbit 
to their mult i-epoch near-infrared RVsH Combining these RVs with the parallax and proper motion 
measured bv lGatewood fc CobanI (|2009l . assuming 4 mas and 4° errors in their proper motion ampli- 
tude and PA, respectively), we derive heliocentric velocities for LP 349-25 of (U, V, W) = (+4.8ib0.8, 
-17.2 ± 1.4, -h5.2 ± 1.3) km s"^ and (+1.9 ± 0.4, -11.2 ± 0.4, -0.5 ± 0.4) km s'^, respectively. 
We adopted the sign convention for U that is positive toward the Galactic center and accounted 
for the errors in the parallax, proper motion, and RV in a Monte Carlo fashion. These space 
motions place LP 349-25 only 0.7a or 0.9a away from the mean of the ellipsoid defined by all 
ot her ultracool dwarfs with UVW measurements. (This sample is described in detail in Section 3.4 
of iDupuv et al.ll2009d .) Using the same method as in our previous work, we have also assessed 
LP 349-25 's m embership in the Galactic thin and thick disk populations using the Besangon model 
of t he Galaxv (IRobin et al.ll2003l ) , finding at most a 0.3% probability of thic k disk membership using 
the lReiners fc Basril (12003) RV and < 0.1% for the IKonopackv et al.l (|201ol l RV. Thus, although the 
derived space motions are different, the conclusion that LP 349-25 is a normal thin disk member is 
unaffected. As a likely thin disk member, the age of LP 349-25 is essentially unconstrained by its 
kinematics. Fina lly, the fact that LP 349-25 is chromospherically active (log(LHa/iboi) = —4.52; 
Gizis et al.ll2000l ) could potentially provide an age constraint, as the activity of M dwarfs change s 
with age. However, the long lifetime of activity in such late-M dwarfs found by IWest et al.l (|2008l ) 
places only a weak constraint of < 9 Gyr for the age of the LP 349-25AB system. Given the mea- 
sured mass and luminosities of LP 349-25AB, it is in fact likely to be much younger, as discussed 
in detail in the following sections. 



For LHS 1901, there is no published radial velocity, precluding the kinematic analysis described 



®The discrepancy between these two RV measurements may be due to the fact that iReiners fc Basril ( 2003 ) 
measured the RV of the primary component, which would ha ve been near the n adir i n its orbital RV curve. Only 
65 days prior to the HIRES measurement (Aphase = 0.02)jjKonopackv et al. (20ldf) measured a primary RV of 
— 11 ± 3 km s~^ , which is more consistent (at 1.6cr) with the IReiners fc Basri (:2009il value. 
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above. The proper motion and parallax of iLepine et al.l (|2009l ) provide a measurement of its 
tangential velocity (T4an = 41.1 it 1.6 ki n s~^), indicating th at LHS 1901 is typical of late-M dwarfs 
within 20 pc (Vtan = 29 it 21 km s~^: Fahertv et al. 20091 ) and not likely to be a member of the 
thick disk or halo populations. Both iReid et al.l ( 20031 ) and iLepine et al.J ( 2009 ) have show n that 
LHS 1901 is not active, which is uncommon for nearby, late-M dwarfs. IWest et al.l (|2008l ) found 
that activity lifetime increases monotonically with M dwarf spectral type, and for M6 and M7 
dwarfs these lifetimes are 7.0 it 0.5 Gyr and S.OI^^q Gyr, respectively. Thus, the lack of activity 
seen in LHS 1901 (M6.5) indicates that its age is likely to be at least as old as 6 Gyr, the 2a lower 
limit for the activity lifetime of M6 and M7 dwarfs. 

The age of th e M2.5 primary star in the G l 569 triple system has been discussed extensively 
in the literature. IZapatero Osorio et al.l (120041 ) find that the space motion of Gl 569A places it 



close to, but not likely a member of, the Ursa Major moving group (300 Myr; ISoderblom et al. 
19931 ). How ever, they suggest th at it may be a member of a stream of young A and F stars 



identified bv IChereul et al.l (119991 ). which have ages of 300-800 Myr as derived from Stromgren 
photometry. ISimon et al.l (120061 ) prefer a younger age of 100-12 5 Myr because Gl 569 A seems to 
lie on the Pleiades sequence in the color-magnitude diagram of iLuhman et al.l (J2005l ). We have 
considered whether the X-ray emission from Gl 569A may also place an independent constraint 
on the a ge of the system. We computed the X-ray flux of Gl 569A from its ROSAT count rate, 
using the lSchmitt et al.l ( 19951 ) conversion factor, finding Fx = (3.03 it 0.30) x 10~^^ erg cm~^ s~^ 
(log Lx = 28.54). This lies between the m edian values of Ly for M dwar fs in the Pleiades (125 Myr) 
and Hyades (625 Myr) as determined by iPreibisch &: FeigelsonI (|2005l ). The X-ray luminosity of 
Gl 569A is higher than ~75% of Hyades objects and is fainter than ~70% of Pleiades objects. This 
suggests an intermediate age of ~300-450 Myr for the Gl 569 system, although its X-ray luminosity 
is consistent with the full range of ~125-625 Myr at la. 



4. Tests of Models 



Directly measured dynamical masses enable unique tests of theoretical models of very low- 
mass stars and brown dwarfs. We consider two independ ent sets of evolutiona ry models from 



the Tucson (JBurrows et al.l 119971 ) and Lyon gro ups f Dusty: IChabrier et al.ll2000ll. F ollowing the 
approach developed in our previous work (e.g., iLiu et al.l l2008l : iDupuy et al.ll2009bl ). we use the 
most directly measured properties of our sample binaries (i.e., total mass and individual component 
luminosities) to derive additional properties (e.g.. Teg) from evolutionary models. These model- 
derived properties are then compared to other available constraints (e.g., from atmospheric model 
fitting). We note that currently available evolutionary models do not incorporate many of the 
theoretical advances used in the latest model atmospheres, such as detailed dust cloud models and 
more accurate molecular and dust opacities. However, our line of analysis is necessarily centered 
on evolutionary models because we directly measure mass and -Lbob not the model atmosphere 
parameters Tefj and log{g). In addition to this analysis, we also explore the approach recently 
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utilized by iKonopackv et alj (J2010l ) in order to directly compare our results to their findings for 
the same models. 



4.1. Model-derived Ages 



As described in detail by lLiu et alJ ()2008l ) and lPupuy et alJ (J2009bl ). the total mass of a binary 
along with its individual component luminosities can be used to estimate the age of the system 
from evolutionary models. This age estimate can be quite precise (!=alO%) when both components 
are substellar, since the luminosities of brown dwarfs depend very sensitively on age. This is the 
case for LP 349-25 AB: the Lyon and Tucson models give consistent ages of 0.127_Q'Q]^y Gyr and 
0.141^gQ]^g Gyr, respectively (Figure[5l Table[7]). The model-derived age of the system is consistent 
with its kinematics (Section 13. 5p . and the implications of such a young age for the LP 349-25AB 
system are addressed in detail in Section [5l 

Because LHS 1901AB is composed of two very low-mass stars whose luminosities remain essen- 
tially constant over their main-sequence lifetime, the model-derived age for this system is not very 

_l_Q '7'7 _l_Q fi'i 

precise. Lyon and Tucson models give consistent ages of 0.28_QQg Gyr and 0.37_q'^5 Gyr, respec- 
tively (Table [8]) . The large upper limits to the age distributions are set by the oldest age for which 
the model grids are computed (10 Gyr). These broad age distributions are consistent with the age 
constraint that comes from the lack of chromospheric activity in LHS 1901 (> 6 Gyr; Section [3. 5p . 
Because the median model-derived ages are quite young, they would be inconsistent with an age as 
old as 6 Gyr, but a more precise parallax measurement is needed to test this. For example, if we fix 
the parallax at its presently measured value but improve its precision from 3.0 mas to 1.0 mas, an 
attainable goal with seeing-limited astrometry, this age discrepancy would rise to > 2a significance 
for the Lyon models. 

Gl 569Bab has a higher mass than LP 349-25AB but lower component luminosities, and thus 
the Lyon and Tucson models give older ages of 0.461q'q7 Gyr and 0.5llQ'gg Gyr, respectively 
(Figure [5l Table [9]). These are consistent with the age constraints fr om Gl 569A discu ssed in 
Section [321 with the exception of the age of 100-125 Myr proposed by ISimon et al.l (J2006l ) based 
on its position on the color-magnitude diagram. Given the unavailability of such data for M dwarfs 
at ~500 Myr, we suggest that Gl 569A may well be consistent with such an older age. Our model- 
derive d age for the Gl 569Bab sys t em is also in agreement with the age estimates from lLane et al.l 
(|200lh and lZapatero Osorio et aD hooj }. 



4.2. Individual Masses 



All of our target binaries have flux ratios near unity, thus enabling robust individual mass 
estimates with only a very weak dependence on model assumptions. For LHS 1901AB, the Lyon 
and Tucson models give mass ratios {q = Mb/Ma) of 0.9581q'q^4 and 0.966^'^Q'Q^g, respectively. Thus, 
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the mo del- derived individual masses (listed in Table [8]) are nearly identical, and their precision is 
dominated by the 1Lii% uncertainty in Mtot- 

Because the model-derived ages for LP 349-25AB are quite young, the inferred mass ra- 
tios are somewhat further from unity than they would otherwise be for the measured 0.133 it 
0.019 dex luminosity ratio. The Lyon and Tucson models give consistent mass ratios of 0.872^'^Q'Q^g 
and 0.863+J];[!}^, respectively. This results in individual masses of 0.064+°;°°| Mq and 0.056 ± 
0.004 Mq for the components of LP 349-25AB (Lyon), placing both below the substellar boundary 



(~ 0.070 M,7^: IChabrier et al.ll2000ll and LP 3 49-25B hkely below the lithium depletion boundary 



0.055-0.065 Mg uChabrier &: Baraffelll997l ). iKonopackv et al.l (|20ld ) reported resolved RV mea- 



surements for LP 349-25AB and, combined with their astrometric orbit, found q = 2.0~^qy. This is 
nominally inconsistent with our model-derived mass ratio at 1.6cr, with their value corresponding 
to the secondary being more massive than the primary. In Figure [U we show their measurements 
along with the RV curve predicted for mass ratios of 0.87 and 2.0, assuming their best-fit center-of- 
mass velocity of —8.0 it 0.5 km s~^. Fixing the orbital elements to the values from the astrometric 
orbit, we find x^ values for the RV measurements that are unreasonably small for both values of 
Q ix^ = 0-5 for q = 2.0 and x^ = 1-7 for q = 0.87; 6 DOF). This implies that the errors on the 
RV measurements are overestimated, and within these large errors both values of the mass ratio 
(0.87 and 2.0) are acceptable. We discuss the plausibility of these different mass ratios in detail in 
Section [5l 

For Gl 569Bab, the Lyon and Tucson models give consistent mass ratios of 0.8661q'q^4 and 
0.886 ln'ni7, respectively. Thi s agrees with the spectroscopically determined mass ratio of 0.711q\3 
from IKonopackv et al.l (|20ld ). Our slightly higher model-derived mass ratio results in somewhat 
more similar masses (0.075ib0.004 and 0.065lQQg4 Mq, Lyon) than implied by the directly measured 
mass ratio (0.082+|];|][]^ and 0.059lo:oo7 Mq). In both cases however, Gl 569Bb is expected to he 
near the mass-limit for lithium burning, which can be tested directly in the future with resolved 
optical spectroscopy of the Li I doublet at 6708 A. In addition, absolute astrometric monitoring of 
the 2.4-yr orbital period binary can readily yield a refined measurement of the mass ratio: 1 mas 
astrometry would give a 5% mass ratio error, ~4x better than spectroscopy. 



4.3. Temperatures and Surface Gravities 

Effective temperature (Tgfj) is one of the most difficult properties to measure directly, as it 
requires a direct measurement of the radius and luminosity. To date, radius measurements remain 
elusive for very low-mass stars and brown dwarfs in the field. Thus, when testing model predictions 
we are restricted to consistency checks between different methods of estimating Tgff . In Section [3. 4^ 
we estimated the temperatures of the primary components of our sample binaries by fitting their 
integrated-light spectra with atmospheric models. The best-fit Drift-PHOENIX and Gaia-Dusty 
model spectra are shown in Figure [7] plotted over the observed spectra. Evolutionary models provide 
a nearly independent temperature estimate by providing a prediction of the radius, which yields 
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Tgff when combined with the measured LboilII 



For LHS 1901A and LHS 1901B, Lyon evolutionary models give effective temperatures of 

50 
40 



2860 =b 50 K and 2820I4Q K, respectively, and the Tucson evolutionary models give 100 K warmer 



temperatures of 2960 ± 30 K and 2930]'^4o ^- Both sets of models give consistent surface gravities 
of log(gf)= 5.2 (cgs) for both components. All NIR atmospheric model fits gave systematically 
higher temperatures of 3000-3100 K and surface gravities near the predicted value (5.0 =b 0.5 dex). 
Given the nominal uncertainty (i.e., the model atmosphere grid step) of 100 K in the best-fit 
Teff values, the largest discrepancy of ~150 K with the Lyon evolutionary models is marginally 
significant. To provide an additional point of comparison f or the Tpff est i mates , we compiled all 



M6.5 and M7 dwarfs with effective tempe ratures derived bvlGautier et al.l (|2007l ) using the nearly 



model-independent infrared fiux method ( Blackwell Sz Shallid 119771 ). We found a mean and stan- 



dard deviation of 2710 ±30 K for these objects, which is lower than Tgff estimates from both classes 
of models by 140-290 K, agreeing somewhat better with evolutionary models. 

For LP 349-25A and LP 349-25B, Lyon models give effective temperatures of 2660 it 30 K and 
2520 lb 30 K, respectively, and the Tucson models give 120 K warmer temperatures of 2780 it 30 K 
and 2640 it 30 K. Lyon models predict surface gravities of about 4.9 (cgs) for both components, 
and Tucson models predict 5.0 (cgs). In comparison, atmospheric model fits of LP 349-25A gave 
2800-3000 K and surface gravities near the predicted value (5.0 ± 0.5 dex)0 Again, the Tucson 
evolutionary model Teff is consistent with that derived from atmospheric models, while the Lyon 
evolutionary models are lower by ~150 K. 

For Gl 569Ba and Gl 569Bb, Lyon models give effective temperatures of 2430 it 30 K and 
2210 lb 30 K, respectively, and the Tucson models give 100 K warmer temperatures of 2530 it 30 K 
and 2300 it 30 K. Lyon models predict surface gravities of about 5.2 (cgs) for both components, 
and Tucson models predict 5.3 (cgs). In comparison, atmospheric model fits of Gl 569Ba gave 
2700-2900 K, with a wide range of surface gravities (log(g) = 4.5 to 6.0). These Teff values are 
much higher than evolutionary models by 150-250 K, and the broad range of surface gravities is in 
reasonable agreement. 

In order to understand why temperature estimates from the two classes of models typically do 
not agree, we show our observed spectra in Figure [8] plotted with Drift-PHOENIX and Gaia-Dusty 
model spectra that have Teff and log(g') values closest to the Lyon Dusty evolutionary model- 



^Nearly independent because evolutionary models use theoretical atmospheres as boundary conditions when eval- 
uating hydrostatic equilibrium and the amount of energy being released. This can result in radius deviations of up 
to ~10% when comparing dustless to extremely dusty model atmospheres at lo w temperatures (Teff < 2000 K), but 
the effect is < 1% at the warmer temperatures of late-M dwarfs (see Figure 2 of lChabrier et al.ll200Gl ). 



*We note that all SED fits including the 0.81 to 0.95 /im optical segment gives a surface gravity of 4.0 (cgs) 
for LP 349-25A, which is highly discrepant with the evolutionary model value of 5.0 (cgs). This would require an 
unrealistic radius of ~0.4 Rq given the known mass of LP 349-25A. Such an unrealistic surface gravity is also found 
from the Ames-Dusty models for LHS 1901A. 
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derived values. This enables us to assess what spectral features prevented our fitting procedure 
from selecting these model spectra, which are very similar to the spectra used for the atmospheric 
boundary condition in the Lyon modelsO The NIR SED fits do not seem to be significantly affected 
by broadband colors (e.g., the better fitting models do not necessarily have J — K colors closer 
to what is observed). Rather, the shapes of the spectra, which are sculpted by broad molecular 
absorption bands, seem to determine the best-fit models. 

• The disagreement in shape is most pronounced at H band, where the observed spectrum is 
much flatter than the low-Tefj, low-gravity model spectra that would better correspond to 
evolutionary model-derived propertie s. This could be due partly to missing FeH opacity in 



2003). 



the rn odels between ~1.5 and 1.7 //m (JReid et al.ll200ll : IWallace k, Hinkldl200ll : ICushing et al 



• The observed X-band shape is also flatter (i.e., the band peak at 2.2 /xm is more suppressed) 
than in the low- Teg, low-gravity model spectra. 

• In J band, the fit seems to be driven by the deep H2O absorption band at 1.33 /xm, which is 
too deep in the low- Toff, low-gravity model spectra. 

• In y band, the FeH 0—0 bandhead at 0.99 //m is too strong in the spectra corresponding to 
evolutionary model-derived properties. This could alternatively be interpreted as the overall 
continuum at Y band being too high in the atmospheric models (e.g., missing opacity from 
wings of resonant alkali lines such as K I) . 

In general, we have found that atmospheric models predict systematically higher effective 
temperatures than evolutionary modelscj This pattern is illustrated in Figure [U which shows the 
Teff estimates from each class of models plotted against each other. A systematic shift of ~250 K 
would bring the models into agreement, with either the atmospheric model estimates being too 
warm or the evolutionary model estimates too cool (i.e., model radii too large by 15-20%) lij We 
note that the observed discrepancy cannot simply be due to the fact that our method relies on fitting 



'The Ames-Dusty model atmospheres provided the spectra used for the boundary condition in the Lyon Dusty 
evohitionary models. However, we choose to show the Gaia-Dusty models for comparison as their SEDs are very 
similar but with updated line lists (Table [5]). This enables us to identify remaining problems and not dwell on issues 
that have already been resolved in the latest generation of that model grid. 

^ iRice et al.1 (2010|) have recently fit model atmospheres to high- and medium-resolution J-band spectra of mid- to 



late-M dwarfs. Most of their sample are young M dwarfs (< 10 Myr), but for their three field M7 to M9 dwarfs they 
find higher T^a values than previously published estimates based on measured luminosities and estimated radii for 
objects of similar spectral type. Thus, their findings are consistent with our results. 

^^ Interestingly, the SED fits that include the wavelength segment from 0.81 to 0.95 /^m (which are less preferred 
for reasons described in Section[3]4]) give Tcft values wlOO K more consist ent with ev olutionary models. This indicates 
that optical diagnostics utilizing TiO bandhead strengths (e.g., see Mo hantv et al . 2004) may be needed to provide 
more accurate T^a estimates than can be obtained from near-infrared spectroscopy alone. 
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the integrated-light spectra of binaries. As discussed in Section 13.41 if the secondary component's 
flux impacted the best-fit model it would have the effect of lowering the derived T^s, whereas we 
find temperatures that are too warm. This means that the discrepancy could only be larger than we 
observe. In addition, even a si gnificant nonsolar n ietallicity for our binaries could not explain the 



observed temperature offset, as lBowler et al.l (|2009l ) showed that models of varying metallicity gave 



the same best-fit T^s to within the model grid step of 100 K for the ultracool subdwarf HD 114762B 
(d/sdM9). 

We have considered what systematic errors could contribute to this discrepancy between evolu- 
tionary and atmospheric models. Currently available evolutionary models use much older versions 
of model atmospheres a s their boundary co r iditio n, and this could potentially result in systematic 



errors in their output. ISaumon &: Marleyl (120081 ) calculated evolutionary models both with and 



without the effects of dust clouds in the atmosphere and found at worst 6% differences in the 
resulting radii, which is much lower than what is needed to account for our observed discrepancy. 
We can also estimate how other input assumptions impact evolutionary model output simply by 
comparing the predictions of the Lyon and Tucson models, which use different helium abundances, 
model atmospheres, and dust treatment. These two sets of models predict radii different by 7- 
10% and temperatures different by ~100 K for the objects in our sample - more than a factor of 
two lower than what is needed to fully account for the observed discrepancy. Presently available 
model atmospheres are more advanced than evolutionary models, but as shown in Figure [7] they 
cannot completely match our observed near-infrared spectra. Therefore, the resulting temperature 
values must harbor systematic errors at some level. Thus, we conclude that the model atmospheres 
must be responsible for at least part of the observed 250 K discrepancy. In fact, we suggest that 
the atmospheric models are more likely than the evolutionary models to be the primary source of 
the discrepancy, since roughly the same Tgfj offset is observed over a wide range of masses, ages, 
and activity levels but the same temperature range. This proposition will be readily testable with 
eclipsing binary measurements given the large implied radius difference between models (15-20%). 



4.4. Near-Infrared Colors 

We computed resolved photometry for our sample binaries in the MKO photometric system to 
compare to Lyon model predictions. For Gl 569Bab, we used our integrated-light MKO photometry 
described in Section 12.11 along with the best available flux ratios from Table [2l For LP 349-25 and 
LHS 1901, we first converted their integrated-light 2MASS photometry to the MKO system by 
deriving correction terms from our near-infrared spectra. The resulting component magnitudes for 
all three binaries are listed in Table HI Near-infrared colors on the MKO photometric system were 
derived from the Lyon models in the same fashion as other properties (e.g., Tefj), and the results 
are listed in Tables [71 [HI and [H In Figures [10] and [TTJ we show the photometry of LP 349-25 AB and 
Gl 569Bab compared to Lyon mass tracks computed for the individual masses of both components. 
The observed JHKL' photometry is significantly discrepant with model tracks at all ages, with 
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models typically predicting ~0. 1-0.2 mag bluer JHK colors than observed and ~0. 1-0.2 mag 
redder K — L' colors. For LHS 1901 AB, the JHK colors predicted by models are in much better 
agreement with the observations. The most discrepant color is H — K, which is only 0.07 mag 
(2.3(t) redder than predicted, and this may be due to missing FeH opacity in the H band. (We do 
not show a corresponding figure for LHS 1901AB as its components are main-sequence stars that 
evolve at essentially constant color.) 



4.5. An Alternative Approach: Model-derived Properties Using Lboi and Teg 



Konopackv et al.l (J2010l ) have recently described an approach to test models of ultracool dwarfs 
with dynamical mass measurements that is somewhat different from ours. They use the measured 
individual luminosities along with Tefj estimates for each component to derive any other property 
(e.g., Mtot or age) from evolutionary models. (Their Tes estimates are derived by fitting model 
atmospheres to resolved broadband photometry of the binaries.) This approach is similar to ours 
in that it uses two quantities to determine a third from evolutionary models, which is necessary 
because brown dwarfs occupy a mass-Lboi-age (or Tpff-Lbni-age) re lation instead of the simpler 
mass-Lboi relation for stars on the main sequence. In lLiu et al.l (|2008l ). we considered in detail such 
model tests using different combinations of these fundamental parameters, and we ultimately chose 
to use mass and Lboi to derive other properties as these are the most directly measured quantities. 
However, any discrepancies we have foun d in the models should al so manifest themselves in the Tgg- 
Lboi approach that has been adopted by iKonopackv et al.l (120101'). In this section, we examine the 
results of analyzing our late-M dwarf measurements using the IKonopackv et al.l (|2010l ) approach. 



We first interp olated the Lyon and Tu cson model grids at individual (Tofj, Lboi) points in 
the same fashion as IKonopackv et al.l (120101 ). with our temperatures coming from the Gaia-Dusty 
atmospheric model fitting of the NIR SEP as descr i bed i n Section 13.41 (This is the exact same 
model atmosphere grid as used by iKonopacky et al.ll2010l .) In order to compute uncertainties in 
the model-derived properties using this approach, we used 10^ luminosity and temperature values 
randomly drawn from normal distributions corresponding to the measurement errors. We assumed 
100 K errors in T^s, equal to the model atmosphere grid step. The resulting evolutionar y model- 
derived masses are plotted in Figure [12] against those derived by IKonopackv et al.l (J20ld ) . 



Our masses are significantly higher than theirs by about a factor of ~2 to 4 because our T^q 
estimates are systematically higher while their luminosity measurements are very similar to ours. 
In order for models to match our much higher temperatures at similar luminosities, the masses 
must be larger and the ages s omewhat older. The bes t illustration of this disagreement is the case 
of 2MASS J22 06— 2047A, a. s IKonopackv et al.l (J2010|) used essentially identical luminosities and 
photometry as iDupuv et al.l (|2009al ). and the co mponent fiuxes are nearly identical, removing any 
ambiguity in fitting its integrated-light spectrum. IKonopackv et al.l (|20ld ) found 2350 ± 80 K from 
fitting its resolved broadband photometry, while our best-fit spectrum from the very same grid of 
model atmospheres has a temperature of 2900±100 K. Therefore, despite the other similarities in our 
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inputs, we found a primary mass of O.OOOlQQQg Mq using the Lyon Dusty evolutionary models, while 
they foun d 0.047ln'n]7 M(T). Ne ither of these masses agrees wi th the measured inass (0 -077ln'n] 7 Mq, 



Table 7 of lDupuv et al.ll2009al ): our mass is higher, and the iKonopackv et al.1 (|201Cl ) mass is much 
lo wer. The differenc e s betw een our derived masses for LP 349-25A and Gl 569Ba are even larger, 
Konopackv et alj ( 20ld ) found extremely low masses for these objects (0.02 ± 0.02 Mq and 



as 



0.02 



h0.02 
-0.01 



Mq, respectively) o 

Fundamentally, both the IKonopackv et al.l ( 20ld ) approach and our standard method to test 
models provide similar information: a consistency check on the temperatures from atmospheric 
and evolutionary models. The comparison presented here highlights how the amplitude, and even 
the sign, of discrepancies found between measured and model-predicted masses depend strongly 
on the input temperature estimates. Using the exact same atmospheric model grids, but different 
methods, our two groups have arrived at values of Tp^ different by f«650 K. This is due to the 
fact that the fitting method of IKonopackv et al.l (|2010l ) relies solely on the broadband colors of 



models, whereas our method is dri ven by the shape of specific spectr al features. As shown in 

b||d), the broadband colors 



Figure [To] and in our previous work (JLiu et al. 



200s 



Dupuy et al. 



2009a 



of models consistently disagree with observations of objects of k nown mass over a broad range of 
spectral types. We also note that a significant disadvantage of the IKonopackv et al.l (|2010l ) Lboi-T'eff 
approach is that evolutionary models occupy a very thin strip in this parameter space as opposed 
to L]^oi~M. This resulted in many of our randomly drawn (Lboh T^g) measurements falling outside 
the range predicted by models (from 1.2% to 99.8% in the worst case), and we simply excluded 
these points from the analysis. 



Konopackv et al.l (|2010l ) interpreted their observed discrepancies with theory as an error in the 



evolutionary model cooling curves, i.e., a mass problem. Using their method, we find an opposite 
mass discrepancy because our Tefj estimates are ~650 K higher than theirs, resulting in higher 
masses. We suggest that the strong sensitivity of this model-testing method to the input T^s 
estimates means that any observed discrepancies more readily identify problems in the way T^s is 
determined from atmospheric models (i.e., a temperature problem), rather than a problem with 
evolutionary models. 



5. The Nature of LP 349-25AB 

It is somewhat surprising that LP 349-25AB h as turned out to be a p air of young (~140 Myr) 
brown dwarfs givei i the lack of detec table lithium (JReiners &: Basrill2009l ) and its M8.0 it 0.5 spec- 



tral classification ( Gizis et al.ll200d ). With a total mass of 0.120lg'QQ7 Mq, one or both of the 



components may even lie below the theoretically predicted lithium depletion boundary at ~0.055- 



( 2OI0I ) imply 
very young ages that are inconsistent with the lack of any low-gravity signatures in these objects' spectra. 



We note that the extremely low model-derived masses of R:i20 Mjup reported bv IKonopackv et al 
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0.065 Mq. In fact, even objects massive eno ugh to deple te lithium take time to achieve the in ternal 
temperature necessary to do so (~ 2 x 10^ K: |Pozidll99ll ). For example. IChabrier et al.1 ( 19961 ) show 
that a 0.070 Mq object takes 220 Myr to destroy 99% of its initial lithium, while at 145 Myr only 
> 0.080 Mq stars are similarly depleted. Using the Lyon and Tucson models, we find predicted 



lithium fractions oi0.68toit and 0.75+°'^^ 



52 for LP 349-25A, respectively, and 0.95l[!;'j'^ 



and 0.96 



-0.03 
-0.08 



for LP 349-25B. Thus, both sets of models predict that the primary may be depleted in lithium, 
but that the secondary should retain nearly all of its initial lithium. The presence of lithium in 
the fainter secondary may be masked by the primary flux in the integrated-light spectroscopy of 
the system, so resolved spectroscopy to detect the Li I doublet at 6708 A is needed to test this 
prediction of brown dwarf evolutionary models directly. 

Evolutionary models make very precise predictions for the age of LP 349-25AB from its mea- 
sured total mass and individual luminosities. Taken together, the Lyon and Tucson models predict 
an age of 140 ±30 Myr. This is consistent with the age of 125 Myr derived for the Pleiades by 
Stauffer et al.l (119981 ) using the lithium depletion boundary. Thus, we might expect LP 34 9 -25AB 
to obey roughly the same boundaries in spectral type and color as found by IStauffer et al.l () 19981 ). 
The lithium boundary in the Pleiades lies between spectral types of M6.5 and M7, with objects 
> M7 having detectable lithium. Although the spectral type determinations for the components 
of LP 349-25AB are hampered by the lack of resolved spectroscopy, the measured integrated-light 
type of M8 .0 =fc 0.5 makes it unl ikely that LP 349-25A is earlier than M7, yet it is clearly depleted 
in lithium. IStauffer et al.l (|l998l ) used spectroscopically derived {R — I) colors to define the lithium 
depletion boundary in the Pleiades, with all ob jects bluer t han 2 .20 mag being depleted in lithium. 
From the integrated-light optical spectruni of iReid et al.l (120031 ). we derived {R — I) = 2.50 mag 
using the same m ethod as 



Stauffer et 



al. 



(Il998l )1^3| 



At this color index, all brown dwarfs in the 



Pleiades sample of IStauffer et al.l (|l998l ) show evidence for the presence of lithium. 



LP 349-25A thus stands out as potentially anomalous compared to determination of the lithium 
boundary in the Pleiades. (Note that our comparison is purely empirical in nature and does not 
depend on any theoretical assumptions regarding lithium depletion in brown dwarfs.) This indicates 
that there may be a problem with the model-derived age for LP 349-25AB. At older ages, the 
lithium depletion boundary moves to lower masses (i.e., later spectral types), so an older age for 
LP 349-25AB would bring it into better agreement with the Pleiades. The age derived in our 
analysis is based solely on the model-predicted Lboi evolution, and for its actual age to be older 
either: (1) the measured parallax would have to be smaller (i.e., LP 349-25AB more distant); or 
(2) the models would need to under-predict the luminosities of LP 349-25AB. In the former case, a 
parallax that is 3.5(T smaller (placing the binary at 14.3 ± 0.3 pc) would increase the total mass to 
0.152 Mq, resulting in a model-derived age 190 it 50 Myr and individual components predicted to 



The method used bv lStauffer et al.l ( 19981 ) was cahbrated against M dwarfs in the Gliese catalog ( Gliese fc Jahreiss 



199lh . and the scatter in (R-1) color at a given M subtype is 0.07 to 0.14 mag. Therefore, we estimate the uncertainties 
in spectroscopically derived values of {R — I) to be «0.10 mag. 
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be nearly fully depleted in lithium. In the latter case, the luminosities we measure would indeed be 
high for their mass, not because the objects are very young but rather because mode l Lbd evolution 



is not correct. Such an error has previously been suggested bv lDupuv et al.l (J2009b) to explain the 
higher-than-predicted luminosities of the benchmark brown dwarf binary HD 130948BC (L4+L4). 
An independent measurement of the parallax for LP 349-25AB and resolved optical spectroscopy 
is needed to better assess this possible discrepancy with brown dwarfs in the Pleiades. 

The youth of LP 349-25AB could potentially be refl ected in low-gravity spectral features such 



as have been noted for other field ultrac ool dwarfs (e.g., iMcGovern et al.ll2004l : lAllers et al 



Kirkpatrick et al.l 120081 : iRice et al.l l20ld ) . The integrated- light optical spectrum from 



2007 



Reid et al 



(|2003l ) and our near-infrared spectrum do not show obvious hallmarks of low surface gravity. How- 



Bihain et al. 



ever, n ear-infrared spectra of M7 and MS dwarfs in the Pleiades recently published by 
(|20ld ) also show no distinguishing characteristics relative to field dwarfs of the same spectral type, 
at least at low resolution, R ~ 50, and modest S/N . We have also checked for possible association 
of LP 349-25AB with young moving groups using our derived (C/, V^ W) helioc entric space velocity 
(Secti on 13.51). We find no suc h association with any of the groups listed by IZuckerman k. Song 
(|2004l ) and lTorreset al.l jiooj)- 



LP 349-25AB is one of the relatively few ultracool dwarfs that d isplays radio emission (~10% 
occurrence rate for >M7 dwarfs: lBergeiil2006l ). iPhan-Bao et al.l (J2007l ) discovered the radio emission 
from LP 349-25AB at 8.5 GHz but with a beam size of 9'.'2 x 8'.'0 were unable to resolve the 



binary to determine whether both components are radio-luminous. lOsten et al.l (|2009l ) presented 
additional unresolved, multi-frequency, multi-epoch observations to better constrain the emission 
mechanism, finding a lack of variability in both the radio fiux and spectral index on both short 
(hours) and long (months) timescales. Their observations spanning 10.7 hours would have easily 
captured multiple rota tions of the emitting component, as the measured t;sin(z) of 56 it 6 km s~^ 
( Reiners &: Basrill20ld ) implies a maximum rotation period of 2.9/ sin(i) hours (assuming an average 
model-derived radius of 0.136 Rq) or a period of 2.6 hours if the rotation axis is co-aligned with 
the binary orbit. LP 349-25AB is the only known radio-luminous ultracool dwarf that d o es no t 
display some level of variability. To explain all observed radio properties, lOsten et al.l (120091 ) 
favored a long-lived (> 0.6 yr), high-latitude (i.e., polar) source emitting gyrosynchrotron radiation. 
Future observations at higher angular resolution using the Extended Very Large Array are needed 
to determine whether only one or both components are radio-luminous. Our dynamical mass 
measurement will be a key input in developing customized models to better understand the origin 
of the magnetic structures generating the unique radio emission of LP 349-25AB. 

Finally, we consider the possibility that LP 349-25AB is actually a higher order multiple sys- 
tem. The resolved colors and magnitudes of the two components are consistent with a simple binary 
as the secondary is slightly fainter and redder than t he primary. As discusse d in Section 14.21 we 
find that the resolved radial velocity measurements of iKonopacky et al.l (|2010l ) are consistent with 
both the model-derived mass ratio of 0.87 and their best fit value of 2.0, in which the secondary has 
twice the mass of the primary. In this scenario, the B component would have to actually comprise 
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two fainter brown dwarfs whose combined luminosity is below that of the single A components 
However, in this case the A component's mass would only be about 0.04 Mq, which is well below 
the lithium depletion limit and thus LP 349-25 should display detectable lithium absorption. In 
addition, dividing the binary into a higher order multiple would necessarily require the age to be 
younger for such lower mass component objects to output the same amount of energy. At younger 
ages, it would be even less likely that the most massive component (whatever its mass) would have 
had time to reach the core temperature needed to destroy its initial lithium. Thus, we conclude 
that given all available constraints LP 349-25AB cannot be a higher order multiple. 



6. Conclusions 

We have determined the orbits of three late-M binaries from Keck NGS and LGS AO orbital 
monitoring. From observations spanning 5.9 years (75%) of the 7.8 yr orbit of LP 349-25AB, we 
have determined a total mass of 0.120lQ'QQy Mq. For LHS 1901AB, we have found a highly eccentric 
(e = 0.830 =b 0.005) 16.1-yr orbit from observations spanning 5.9 years (36% of the orbit, including 
periastron passage), deriving a dynamical mass of 0.194_o'o2i -^0- We have also determined a new 
orbit for Gl 569Bab based solely on astrometrically well-calibrated instruments (Keck/NIRC2 and 
JfST/STIS), finding a larger semimajor axis and thus larger dynamical mass (0.140lQ'QQg Mq) than 
previous work. For all binaries, the most significant contribution to their mass uncertainties are 
the errors in their parallaxes. 

Despite LP 349-25AB's integrated-light spectral type of M8 and lack of detectable lithium 
absorption, this binary has turned out to be pair of brown dwarfs. To match the observed total 
mass and component luminosities, evolutionary models predict that the system must be quite young 
(140 lb 30 Myr) and that at least the secondary component should have retained most of its initial 
lithium. If the model-derived age for LP 349-25AB is correct, it is nearly the same age as the 
Pleiades; however, it appea rs to be discrepant w ith the boundary between l ithium-depleted and 



lithium -bearing objects from Stauffer et al.l (|l998l ). If the parallax measured bv lGatewood &: Coban 



(|2009l ) is accurate, this disagreement with the Pleiades indicates a possible problem with the model- 
derived age, which in our analysis is based solely on the predicted Lboi evolution. An older age for 
LP 349-25AB would remedy this problem, as the lithium depletion boundary moves to later spectral 
types at older ages. This solution would require t hat models uiider-pre dict the luminosities of brown 



dwarfs, which has previously been suggested by iDupuy et al.l (J2009bl ) to explain the higher-than- 



predicted luminosities of the benchmark brown dwarf binary HD 130948BC (L4-I-L4). We also 
find that the lack of lithium depletion in at least the primary component effectively rules out the 
possibility that LP 349-25AB is a higher order multiple. 



^^Another possibility resulting in g = 2.0 that can be rul ed out is that the seco ndary component of LP 349-25AB 



is more luminous than the primary. As shown in Figure 1 of lBurrows et al.l (|200l|), a secondary with a putative mass 



of 0.04 M© would never outshine the primary, even during its deuterium-fusing phase at «2.5 Myr. 
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As in our previous work (e.g., iLiu et al.l 120081 : iDupuv et al.l l2009bl ) . we have used measured 
dynamical masses to test the predictions of evolutionary and atmospheric models. Properties such 
as Tgff and log{g) are derived from evolutionary models using only Mtot and the resolved component 
luminosities. We have also fit model atmospheres to the integrated-light 0.95-2.42 /im spectra as 
an independent meth od for estimating Tgff . Including our previous work on 2MASS J2206— 2047AB 
(jDupuy et al.ll2009al ). there is now a sample of four late-M binaries with precise dynamical masses 
and near-infrared spectra enabling a consistency check between evolutionary model-derived proper- 
ties and model atmospheres. All four model atmosphere grids we tested give effective temperatures 
~250 K warmer than predicted for our objects of known mass from Lyon Dusty evolutionary model 
radii (including the Ames-Dusty model atmospheres which were used as the boundary condition 
for the Dusty evolutionary models). 

Although model atmospheres now incorporate considerably more advanced theory than the 
decade-old evolutionary models, we propose that such a large offset cannot be entirely explained 
by problems with model radii. Instead, model atmospheres are likely to be the major source of 
this discrepancy given that: (1) they must harbor systematic errors at some level, since they do 
not completely match our observed spectra; and (2) roughly the same Tgfj offset is observed over 
a narrow range of Teg but a wide range of masses, ages, and activity levels. Such a discrepancy 
would be difficult to explain in the context of evolutionary models, although they may also harbor 
systematic errors that partially contribute to this inconsistency. Directly measured radii are needed 
to determine which set of models correctly predicts Tefj (if either do). This can be readily tested 
with future discoveries of late-M eclipsing binaries, as a 250 K offset corresponds to a substantial 
radius difference (15-20%). 

In add ition to our standard liii e of analysis, we have also explored the approach to testing mod- 
els used bv lKonopackv et al.l (J2010l ) for the late-M binaries common to our two samples. Our model 
atmosphere fitting of integrated-light spectra yields higher Tefj estimates than their approach of 
fitting resolved broadband photometry. Consequently, we find evolutionary model-derived masses 
from their approach that are a factor of ~2 to 4 higher than theirs, resulting in an opposite discrep- 
ancy with measured dynamical masses. This illustrates that their method sensitively relies upon 
the input T^s, a property that cannot be determined directly without radius measurements. Thus, 
this approach may have limited utility in testing evolutionary models, as observed discrepancies 
are more readily caused by problems with Tpff estimates from mod el atmospheres. To summarize 
how our conclusions contrast with those of lKonopackv et al.l (|2010l ): (1) They found estimates for 
Teff from model atmospheres that are ~400 K lower than our evolutionary model-derived values, 
whereas we find T^s estimates that are ~250 K higher. (2) They interpreted these results in terms of 
mass, concluding that systematic errors in the evolutionary model cooling curves were responsible 
for observed discrepancies. We propose instead that the discrepancy is largely caused by systematic 
errors in the model atmospheres. 

As the sample of ultracool dwarfs with dynamical mass measurements grows, we are beginning 
to probe the diversity of the field population, which comprises both very low-mass stars and brown 
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dwarfs. At the cool temperatures of late-M dwarfs, field objects can span nearly a factor a two in 
mass depending on their age, which makes mass measurements crucial for breaking this degeneracy. 
At cooler temperatures there are even fewer mass measurements presently available, and the effects 
of complex processes such as atmospheric dust formation become more important. Therefore, 
continued orbital monitoring will be imperative to bring similar tests of models as we report here 
into the much cooler regimes of the L and T dwarfs. 
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Fig. 1.— Relative astrometry for LP 349-25AB, LHS 1901AB, and Gl 569Bab along with the best- 
fitting orbits. The dotted line and arrow indicate the time of periastron passage, the dashed line 
shows the line of nodes, and the empty circles show predicted future locations. Error bars are 
typically smaller than the plotting symbols. All orbits are sufficiently well constrained that the 
uncertainties in the total masses are dominated by the parallax errors. 
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Fig. 2.— Measurements of the projected separation and PA of LP 349-25AB (left), LHS 1901AB 
(middle), and Gl 569Bab (right). The best-fit orbits from our MCMC analysis are shown as sohd 
Hnes. The bottom panel of each plot shows the observed minus computed (O — C) measurements 
with observational error bars. 
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Fig. 3. — The covariance between the orbital parameters for LP 349-25AB (top), LHS 1901AB 
(middle), and Gl 569Bab (bottom) shown in grayscale using all 2 x 10^ saved steps from our MCMC 
chains of length 2 x 10^. For LP 349-25AB and LHS 1901 AB the correlation between P and a 
enables the total mass to be determined more precisely than from simple propagation of errors, 
as illustrated by lines drawn demarcating the 3o" range for the total mass (prior accounting for 
the distance uncertainty). The double-peaked nature of the MCMC parameter distributions for 
LHS 1901AB and Gl 569Bab also reveal slightly degenerate orbit solutions. 
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Fig. 4. — Contours of the x^ values from model atmosphere fits to our sample binaries' observed 
integrated-light 0.95-2.42 ^m spectra (drawn at 1.02, 1.1, 1.3, 1.5, 1.7, 2.0, 2.5, and 3.5 times 
the minimum x^)- 2^eff is generally better constrained than log{g), as indicated by contours that 
are elongated in the log{g) direction. White space shows where models are not present in our 
grids. For comparison, the gold stars show Lyon Dusty evolutionary model-derived properties, 
with log{g) in general agreement but Teg systematically lower than the best-fit model atmospheres 
(see Section [43]) . 
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Fig. 6. — Radial ve locities of LP 349-25A (blue) and LP 349-25B (red) as measured by 



Konopacky et al.l (120101 ) shown with curves predicted by our best-fit orbit, assuming two different 
mass ratios {q = M2/M1). Combi ning our measured lumi nosity ratio and total mass, evolutionary 
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Fig. 7. — Integrated-light SpeX SXD spectra (black) shown with the best fitting atmospheric model 
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is labeled indicating its effective temperature, with surface gravity in parentheses. We performed 
separate fits for the full NIR SED (left) and individual bands (right). All spectra are smoothed by 
10 pixels Rr: 0.003 /um for display only. 
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Fig. 8. — Integrated-light SpeX SXD spectra (black) shown with the atmospheric model spec- 
tra corresponding to the values of T^s ci^d log(5') derived from Lyon Dusty evolutionary models. 
Each model spectrum (brown) is labeled indicating its effective temperature, with surface gravity 
in parentheses. All spectra are smoothed by 10 pixels ~ 0.003 ^m for display only. The most 
prominent deviations in the model spectra are: the exaggerated shapes of H and K bands; the 
overpredicted depth of H2O absorption at 1.33 /im; and the strength of the FeH bandhead at 
0.99 fiui. Without directly measured radii we cannot determine if the problem lies in the model 
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typically lie above (or to the left of) the line of equality. An offset of 250 K can largely account for 
the observed discrepancies (dashed line) , implying that either evolutionary model estimates are too 
cool (i.e., radii too large by 15-20%) or that atmospheric model estimates are too warm by 250 K. 
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Fig. 10. — Color-magnitude diagrams showing the measured photometry of LP 349-25A (top) and 
LP 349-25B (bottom) compared to Lyon evolutiona ry tracks (all photom etry on the MKO system). 
The solid lines are isomass tracks from the Dusty (jChabrier et al.l 120001 ) models with dotted lines 
encompassing the la mass uncertainties, and open squares demarcate ages of 0.1, 0.2, and 0.3 Gyr. 
Field dwarfs with parallax measurements are shown for comparison as filled gray circles. LP 349- 
25A and B are somewhat redder (0.1-0.2 mag) than evolutionary model tracks for JHK colors, 
while their K — V colors are 0.1-0.2 mag bluer than models. 
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Fig. 11. — Same as Figure do] but for Gl 569Bab. 
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Fig. 12. — Comparison of the primary component masses derived from evolutionary models using 
the Lboi-^eff approach by iKonopackv et al.l ((20101) and by us. The masses we derive using this 
approach are systematically higher by a factor of ~2 to 4 (left) because our model atmosphere 
fitting of integrated-light spectra yielded ~650 K higher temperatures than their fitting of resolved 
broadband photometry (right). We suggest that the sensitivity of this model testing approach to 
the input T^s estimates makes it more likely to identify problems with model atmospheres than 
evolutionary models. 
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Table 1. Keck AO Observations 



Date 


Time 


Airmass 


Filter 


FWHM'' 


Strehl ratio'^ 


(UT) 


(UT) 






(mas) 




LP 349-25AB 


2008 Jan 16<= 


05:56 


1.280 


Ks 


57 ±3 


0.23 ± 0.05 


2008 Jun 30<^ 


13:20 


1.314 


J 


55 ±4 


0.026 ± 0.009 




13:29 


1.274 


H 


48.6 ±0.8 


0.123 ±0.012 




13:12 


1.350 


Ks 


58 ±4 


0.19 ±0.04 


2008 Aug 20*^ 


14:21 


1.069 


Ks 


50.1 ±0.5 


0.448 ± 0.007 


2008 Sep 9*^ 


13:51 


1.169 


L' 


80.2 ±0.9 


0.78 ±0.03 


2009 Sep 28^^ 


13:17 


1.328 


Ks 


54 ±3 


0.29 ±0.05 


2009 Doc IS'^ 


05:22 


1.479 


K 


55.8 ±0.9 


0.33 ±0.03 


2010 May 22<^ 


14:54 


1.721 


K 


58 ±3 


0.36 ± 0.04 






LHS 


1901AB 






2008 Jan 15*^ 


11:05 


1.128 


Ks 


49.0 ±1.9 


0.411 ±0.012 


2008 Sep 9*^ 


15:13 


1.411 


J 


39.8 ±1.8 


0.087 ±0.015 




14:59 


1.475 


H 


42.6 ±1.3 


0.217 ±0.014 




15:17 


1.393 


Ks 


46.7 ±1.1 


0.47 ±0.03 


2009 Sep 28'' 


15:37 


1.150 


Ks 


51 ±3 


0.30 ±0.07 


2009 Dec le'^ 


08:45 


1.428 


Ks 


53.7 ±1.7 


0.30 ± 0.04 


2010 Mar 22*^ 


05:59 


1.096 


J 


41 ±5 


0.05 ± 0.02 




05:55 


1.094 


H 


40.1 ±0.8 


0.23 ±0.03 




05:52 


1.093 


K 


50.3 ±0.2 


0.43 ±0.02 




06:04 


1.097 


L' 


82.7 ±0.4 


0.63 ± 0.07 


Gl 569Bab 


2004 Dec 24^'"^ 


16:22 


1.270 


H 


46.8 ±0.5 


0.17 ±0.04 


2005 Feb 25''''* 


12:08 


1.298 


Hcout 


39.0 ±1.1 


0.23 ±0.03 


2008 Jan 16^ 


16:33 


1.036 


^cont 


36.9 ±0.8 


0.31 ±0.05 


2009 Apr 29*= 


12:24 


1.089 


"cont 


50.0 ±1.3 


0.39 ±0.07 


2009 May 29*^ 


10:56 


1.159 


-ft'cont 


52.0 ±0.7 


0.58 ± 0.04 


2010 Mar 22*^ 


15:48 


1.235 


K 


48.6 ±1.2 


0.50 ±0.06 


2010 May 23*^ 


11:10 


1.133 


J 


37 ±2 


0.067 ±0.011 




11:16 


1.148 


H 


40.6 ±1.6 


0.19 ±0.06 




11:23 


1.166 


K 


49.8 ±0.6 


0.44 ± 0.04 




11:32 


1.193 


L' 


82.4 ±0.5 


0.76 ± 0.02 



^Strehl ratios and FWHM were computed using the publicly available rou- 
tine NIRC2STREHL. Errors are the rms of individual dithers. 

'^NGS AO observations. 

"^LGS AO observations. 

"^Data originally published bv lSimon et al.l 1120061) and reanalyzed by us. 
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Table 2. Best-Fit Binary Parameters 



Epoch (UT) 


Instrument 


Filter 


p (mas) 


PA(°) 


Am (mag) 


LP 349-25AB 


2004 Jul 3 


CFHT/PUEO'^''' 


K' 


125 ± 10 


12.7 ±2.0 


0.26 ±0.05 


2004 Sep 26 


VLT/NACO''''' 


H 


107 ±10 


7.1 ±0.5 


0.38 ±0.05 


2006 Jul 13 


VLT/NACO^ 


Ks 


105.6 ±0.3 


247.3 ±0.2 


0.33 ±0.03 


2006 Sep 21 


VLT/NACO'' 


Ks 


116.6 ±0.6 


240.4 ±0.2 


0.31 ±0.03 


2006 Oct 9 


VLT/NACO'' 


Ks 


118.9 ±0.6 


238.9 ±0.3 


0.285 ±0.013 


2006 Nov 11 


VLT/NACO'' 


Ks 


123.3 ±0.4 


236.0 ±0.2 


0.315 ±0.019 


2006 Dec 24 


VLT/NACO'' 


Ks 


128.2 ±0.7 


232.7 ±0.3 


0.276 ±0.018 


2008 Jan 16 


Keck/NIRC2'' 


Ks 


137.3 ±0.4 


207.92 ±0.09 


0.315 ±0.011 


2008 Jun 30 


Keck/NIRC2 


J 


115.2 ±0.3 


194.4 ±0.4 


0.35 ± 0.03 




Keck/NIRC2'' 


H 


114.87 ±0.12 


194.71 ±0.13 


0.326 ±0.011 




Keck/NIRC2 


Ks 


115.2 ±0.5 


194.64 ±0.13 


0.318 ±0.007 


2008 Aug 20 


Keck/NIRC2'' 


Ks 


105.9 ±0.2 


189.71 ±0.11 


0.314 ±0.007 


2008 Sep 9 


Keck/NIRC2'' 


V 


102.45 ± 0.19 


187.46 ±0.10 


0.222 ±0.005 


2009 Sep 28 


Keck/NIRC2'' 


Ks 


71.1 ±0.3 


98.3 ±0.6 


0.24 ±0.03 


2009 Dec 15 


Keck/NIRC2'» 


K 


83.4 ±0.3 


81.6 ±0.4 


0.38 ±0.06 


2010 May 22 


Keck/NIRC2'» 


K 


112.6 ±0.4 


59.76 ±0.13 


0.307 ± 0.008 






LHS 


1901AB 






2004 Jan 8 


CFHT/PUEO''''^ 


K' 


275 ±5 


208.0 ±0.5 


0.13 ±0.03 


2005 Apr 27 


CFHT/PUEO'^'^ 


K' 


204 ±5 


215.0 ±0.5 


0.07 ± 0.03 


2005 Oct 14 


CFHT/PUEO'''= 


H 


174 ±5 


219.6 ±0.5 


0.14 ±0.05 


2008 Jan 15 


Keck/NIRC2^ 


Ks 


60.5 ±1.5 


308.6 ±0.8 


0.16 ±0.06 


2008 Sep 9 


Keck/NIRC2'» 


J 


57.4 ± 0.6 


1.0 ±0.8 


0.31 ±0.16 




Keck/NIRC2 


H 


57.2 ±1.3 


1.2 ±1.8 


0.28 ±0.13 




Keck/NIRC2 


Ks 


57.8 ±1.0 


1.5 ±0.8 


0.15 ±0.04 


2009 Sep 28 


Keck/NIRC2'' 


Ks 


177.90 ±0.14 


179.54 ±0.03 


0.098 ± 0.007 


2009 Dec 16 


Keck/NIRC2'' 


Ks 


206.78 ±0.13 


181.31 ±0.03 


0.094 ± 0.003 


2010 Mar 22 


Keck/NIRC2 


J 


236.9 ±0.4 


183.12 ±0.03 


0.111 ±0.016 




Keck/NIRC2 


H 


236.9 ±0.4 


183.11 ±0.04 


0.115 ±0.009 




Keck/NIRC2'' 


K 


237.0 ± 0.3 


183.10 ±0.04 


0.107 ±0.007 




Keck/NIRC2 


L' 


237.3 ±0.7 


183.15 ±0.05 


0.099 ± 0.003 


Gl 569Bab 


2002 Jun 26 


HST/STIS'' 


F28X50LP 


97.3 ±1.3 


94.0 ±1.3 


0.99 ±0.03 


2004 Dec 24 


Keck/NIRC2=' 


H 


93.9 ±0.3 


110.8 ±0.4 


0.54 ± 0.02 


2005 Feb 25 


Keck/NIRC2'' 


Hcoat 


85.2 ±0.8 


132.9 ±0.7 


0.55 ±0.08 


2008 Jan 16 


Keck/NIRC2'' 


.^cont 


61.8 ±0.7 


272.5 ±1.4 


0.58 ±0.16 


2009 Apr 29 


Keck/NIRC2'' 


K'cont 


100.2 ± 0.6 


66.5 ±0.6 


0.55 ±0.07 


2009 May 29 


Keck/NIRC2'' 


-K'cont 


100.5 ±1.8 


75.0 ±1.7 


0.59 ±0.07 


2010 Mar 22 


Keck/NIRC2'' 


K 


55.8 ±0.3 


206.5 ±1.1 


0.473 ±0.010 


2010 May 23 


Keck/NIRC2 


J 


58.5 ±1.1 


265.9 ±1.0 


0.45 ±0.13 




Keck/NIRC2 


H 


58.8 ±0.5 


267.6 ±0.6 


0.46 ± 0.06 




Keck/NIRC2'» 


K 


59.9 ±0.4 


268.1 ±0.4 


0.49 ± 0.04 




Keck/NIRC2 


L' 


59.9 ±0.6 


267.6 ±0.5 


0.49 ±0.03 


^■Used in the orbit fit. 










''Measurements from iForveille et alj 1120051"). 
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'^Measurements from lMontagnier et alj l l200(i) . 
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Table 3. Derived Orbital Parameters^ 





MCMC 




MPFIT 


Parameter 


Median 68.3% c.l. 


95.4% c.l. 


routine 


LP 349-25AB 



Seniimajor axis a (mas) 
Orbital period P (yr) 
Eccentricity e 
Inclination i (°) 

Time of periastron passage Tq ~ 
PA of the ascending node Q {°) 
Argument of periastron ui {°) 
Total mass (M©); fitted<: 
Total mass (Mq): final"* 
X^ (21 degrees of freedom) 



2454860.5*^ (JD) 



146.7 

7.76 

0.051 

117.24 



35.95 

250 

0.1205 

0.120 



-0.6, 0.6 

-0.04, 0.04 

-0.003, 0.003 

-0.14, 0.14 

-24, 26 

-0.12, 0.12 

-3, 4 

-0.0007, 0.0007 

-0.007, 0.008 



-1.1, 1.2 

-0.07, 0.08 

-0.005, 0.006 

-0.3, 0.3 

-50, 50 

-0.23, 0.23 

-6, 7 

-0.0013, 0.0013 

-0.014, 0.017 



146.7 ±0.6 
7.76 ±0.04 
0.051 ±0.003 
117.23 ±0.13 
-2 ±25 
35.95 ±0.11 
250 ±3 
0.1205 ±0.0006 
0.120 ±0.008 
18.1 



LHS 1901AB 



Semimajor axis a (mas) 288 —5, 5 —8, 8 288 ± 3 

Orbital period P (yr) 16.1 -0.5,0.5 -0.8,0.8 16.1 ±0.3 

Eccentricity e 0.830 -0.005, 0.005 -0.008, 0.008 0.830 ± 0.003 

Inclination^' j (°) 72.1 -0.3,0.3 -0.5,0.5 72.12 ±0.18 

Time of periastron passage To - 2454765.0== (JD) 0.0 -2.5,2.5 -4,4 0.0 ±1.6 

PA of the ascending node Q (°) 182.0 -0.2, 0.2 -0.4, 0.4 181.99 ± 0.15 

Argument of periastron LJ (°) 224.8 -0.6,0.6 -1.1,1.1 224.8 ± 0.4 

Total mass (Mq): fitted<= 0.1944 -0.0028,0.0028 -0.005,0.005 0.1944 ± 0.0018 

Total mass (Mq); finald 0.194 -0.021,0.025 -0.039,0.053 0.194 ±0.023 

X^ (9 degrees of freedom) ■ ■ ■ ■ ■ ■ ■ ■ • 6.8 



Gl 569Bab 



Semimajor axis a (mas) 95.6 -1.0,1.1 -1.8,2.0 95.7 ±0.7 

Orbital period P (yr) 2.367 -0.003, 0.003 -0.006, 0.006 2.367 ± 0.002 

Eccentricity e 0.316 -0.005, 0.005 -0.009, 0.009 0.316 ± 0.003 

Inclination^' j (°) 35.0 -1.1,1.1 -2.2,2.0 35.1 ±0.7 

Time of periastron passage To - 2455290.5^ (JD) 0.0 -1.5,1.5 -2.5,2.5 0.0 ± 0.9 

PA of the ascending node Q (°) 324.8 -2.0,2.0 -3.3,3.3 324.9 ± 1.2 

Argument of periastron LJ (°) 257.9 -2.0,2.0 -3.0,3.3 257.8 ± 1.2 

Total mass (Mq): fitted<: 0.140 -0.004,0.005 -0.007,0.009 0.1408 ± 0.0028 

Total mass (Mq): final"* 0.140 -0.008,0.009 -0.015,0.018 0.141 ±0.008 

X^ (9 degrees of freedom) ■ ■ ■ ■ ■ ■ ■ ■ ■ 6.3 



^Note that the uncertainties quoted here are "single parameter" errors (i.e., Ax^ = 1) and thus are only valid when a 
single parameter is of interest. 

''2009 January 27 00:00:00.0 UT 

"^The "fitted" total mass represents the direct results from fitting the observed orbital motion without accounting for 
the parallax error. 

''The "final" total mass includes the additional error in the mass due to the error in the parallax. 

==2008 October 25 12:00:00.0 UT 

<=2010 April 4 00:00:00.0 UT 



Table 4. Measured Properties of Target Binaries^ 



LP 349-25 



LHS 1901 



Gl 569B 



Property 



Ref. 



Rcf. 



Ba 



Bb 



Ref. 



A/tot (Mq) 


0.120 


+0.008 
-0.007 


Semimajor axis (AU) 


1.94 ± 0.04 


d (pc) 


13.2 ±0.3 


Spectral type 


M7.5±1.0 


M8.0ti'o 


J (mag) 


11.154 ±0.025 


11.506 ±0.028 


H (mag) 


10.606 ±0.023 


10.932 ± 0.023 


K (mag) 


10.151 ±0.018 


10.458 ± 0.018 


V (mag) 


9.80 ±0.07 


10.02 ± 0.07 


J-K (mag) 


1.00 ±0.03 


1.05 ±0.03 


H-K (mag) 


0.46 ±0.03 


0.47 ±0.03 


J-H (mag) 


0.55 ±0.03 


0.57 ±0.04 


K - L' (mag) 


0.35 ±0.07 


0.44 ±0.07 


Afj (mag) 


10.55 ±0.05 


10.90 ±0.05 


Mh (mag) 


10.01 ±0.05 


10.33 ± 0.05 


Mif (mag) 


9.55 ±0.05 


9.86 ±0.05 


Mu (mag) 


9.20 ±0.08 


9.42 ±0.08 


log(ibol/-f'0) 


-3.041 ±0.024 


-3.175 ± 0.026 


Alog(Lboi) 


0.133 ±0.019 



1,2 

1,2 

2 

1 

1,7 

1,7 

1,7 

1,8 

1,7 

1,7 

1,7 

1,7,8 

1,2,7 

1,2,7 

1,2,7 

1,2,8 

1 

1 



n 104+0.025 

u.iy4_Q 021 
3.70 ±0.16 
12.9 ±0.5 



M6.5tJ;5 
10.631 ±0.020 
10.197 ±0.016 
9.796 ±0.019 

0.83 ±0.03 
0.40 ±0.02 
0.43 ±0.03 

10.09 ±0.09 
9.65 ±0.09 
9.25 ±0.09 



M6.5tJ;° 
10.743 ±0.020 
10.312 ±0.016 
9.904 ±0.019 

0.84 ±0.03 
0.41 ±0.03 
0.43 ± 0.03 

10.20 ±0.09 
9.77 ±0.09 
9.36 ±0.09 



-2.95 ±0.04 -2.99 ±0.04 

0.044 ±0.015 



1,3 
1,3 

3 

1 

1,7 
1,7 
1,7 

1,7 
1,7 
1,7 

1,3,7 
1,3,7 
1,3,7 

1 
1 



n 140+0.009 

U.i4U_Q (JQg 

0.923 ± 0.018 
9.65 ±0.16 



M8.5±0.5 

11.27 ±0.06 

10.67 ±0.04 

10.16 ±0.03 

9.46 ±0.10 

1.11 ±0.07 

0.50 ±0.05 

0.61 ±0.07 

0.70 ±0.10 

11.35 ±0.07 

10.74 ±0.06 

10.24 ±0.05 

9.54 ±0.10 

-3.424 ±0.019 



M9.0±0.5 

11.78 ±0.06 

11.21 ±0.04 

10.64 ±0.03 

9.95 ±0.10 

1.15 ±0.07 

0.57 ±0.05 

0.58 ± 0.07 

0.68 ±0.11 

11.86 ±0.07 

11.28 ±0.05 

10.71 ±0.05 

10.03 ±0.11 

-3.623 ± 0.020 



0.199 ±0.016 



1,4 
1,4 

4 
5,6 

1 

1 

1 
1,9 

1 

1 

1 

1,9 
1,4 
1,4 I 

' 00 
1,4,^ 

1 

1 



^AU near-infrared photometry on the MKO system. 



Refe renc es. — (1) This wo rk; (2 llGatewood fc CobanI ||2009|' ): (3')lLepine et a.1.1 1|2009|) : (4) Ivan LeeuweiJ ll2007l l: (51 iHenrv 
l l2004h : (7) ICutriet ai]||2003l '): (81 IColimowski et al.1 l l2004h : ('9'l lForrest et al.l lll988l '). 



KirkpatricM lll99Cl 1: (61 IZapatero Osorio et al.l 
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Table 5. Description of Model Atmosphere Grids 



Grid name 


Grid 




Ref. 


PHOENIX-Gaia 


BH05 


Ames-Dusty 


AOl 


Gaia-Dusty 


BIO 


Drift-PHOENIX 


W09 



Molecular line lists 



H2O 



TiO 



FeH 



CrH 



Element 
abundance 



Dust treatment 



Ames Ames DA99 F99 GN93 

Ames Ames PD93 F99 GN93 

BT06 Ames DOS B02 A05 

Ames Ames B06 B06 GN93 



dust disappears after forming ("Cond") 
dust stays where it forms ("Dusty") 
dust stays where it forms ("Dusty") 
noncquilibrium cloud model ("Drift") 



Ref erences. — Ames: IPartridge fc Schwenk3 lll997h. ISchwenkj lll998h : AOl: lAUard et all ll200lh : A05: lAsplund et~all 
ll2005ll: B02: iBurrows et alj l|2002h: B06: iBernathT l|2006l'); BIO: T Barm an (2010, private communication); BH05: 
iBrott fc Hauschildj l|2005h : BT06: iBarber et alj ll2006h : D03: iDulick et all ll2003h: DA99: Da.vis fc A Uard (1 999. private com- 
munic ation); F99: R. D. Freedrn an (1999, private communication); GN93: iGrevesse fc Noelj l ll993l) : PD93: IPhillips fc David 
l|l993l l: W09: IWitte et alJ l|200gl ). 
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Table 6. Best-fit Atmospheric Models 





PHOENIX-Gaia 


Ames 


-Dusty 


Gaia- 


■Dusty 


Drift-PHOENIX 


Spectral range 


Toff (K) 


log(g) 


Tcff (K) 


log(5) 


Teff (K) 


log(5) 


Tcff (K) 


log(ff) 


LP 349-25 


NIR (0.95-2.42 fim)'' 


2800 


4.5 


2800 


5.5 


2800 


5.0 


2800 


5.5 


All (0.81-2.42 ^m) 


2600 


4.0 


2600 


4.0'= 


2600 


4.0 


2600 


4.0 


Y (0.95-1.12 ^m) 


2800 


4.5 


2800 


5.5 


2800 


4.5 


2700 


5.0 


J (1.10-1.34 ^m) 


2900 


5.5'= 


2900 


5.5 


2900 


5.5 


2800 


5.5 


H (1.40-1.80 ^m) 


2900 


5.5'= 


2900 


5.5 


2900 


5.5 


2800 


6.0'= 


K (1.90-2.40 ^m) 


3000 


5.5'= 


3000 


5.5 


2900 


5.5 


2900 


6.0*= 


LHS 1901 


NIR (0.95-2.42 /iin)=' 


3000 


5.0 


3000 


5.0 


3000 


5.0 


3000 


5.0 


All (0.81-2.42 fim) 


2900 


4.5 


2900 


4.0'= 


2900 


4.5 


2900 


5.0 


Y (0.95-1.12 ^m) 


3100 


4.5 


3000 


5.0 


3000 


4.5 


3000 


5.0 


J (1.10-1.34 ^m) 


3000 


5.5'= 


3000 


5.5 


3000 


5.5 


3000 


5.5 


H (1.40-1.80 ^m) 


3000 


5.5'= 


3000 


5.5 


3000 


5.5 


3000 


6.0'= 


K (1.90-2.40 urn) 


3100 


5.5'= 


3100 


5.5 


3100 


5.5 


3000 


5.5 


Gl 569B 


NIR (0.95-2.42 /iin)=' 


2800 


4.5 


2800 


5.5 


2800 


5.0 


2800 


5.5" 


All (0.81-2.42 ^m) 


2700 


4.5 


2800 


5.5 


2700 


5.0 


2700 


5.0" 


Y (0.95-1.12 ^m) 


2900 


5.0 


2900 


6.0'= 


2800 


5.0 


2700 


5.5 


J (1.10-1.34 ^m) 


2800 


5.5'= 


2700 


5.5 


2700 


5.5 


2700 


5.5 


H (1.40-1.80 ^m) 


2900 


5.5'= 


2800 


6.0*= 


2800 


6.0*= 


2800 


6.0*= 


K (1.90-2.40 pm) 


2900 


5.5'= 


2900 


5.5 


2800 


5.5 


2800 


6.0'= 


2MASS J2206-2047 


NIR (0.95-2.42 /iin)=' 


2900 


5.0" 


2900 


6.0'= 


2900 


5.5 


2900 


5.5 


All (0.81-2.42 ^m) 


2800 


4.5 


2800 


6.0'= 


2700 


5.0 


2800 


5.5 


Y (0.95-1.12 ^m) 


2900 


5.0 


2900 


6.0*= 


2800 


5.0 


2700 


5.5 


J (1.10-1.34 ^m) 


2900 


5.5'= 


2800 


5.5 


2800 


5.5 


3000 


6.0'= 


H (1.40-1.80 ^m) 


2900 


5.5'= 


2900 


6.0*= 


2900 


5.5 


2900 


6.0*= 


K (1.90-2.40 ^m) 


3000 


5.5'= 


3000 


5.5 


2900 


5.5 


2900 


6.0*= 



''The NIR results are our preferred choice for model atmosphere fitting, as discussed in Section 13.41 

'=The best-fit value is at the edge of the model grid. 

"In these rare cases (~3% of all the tabulated results), a significant x^ minimum developed that was very discrepant from 
the prevailing best-fit models at «2700-2900 K. Thus, we selected the second best fitting model rather than: 2000 K/6.0 for 
Gl 569B Drift-PHOENIX; 3300 K/3.5 for 2MASS J2206-2047 PHOENIX-Gaia. 
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Table 7. Evolutionary Model-derived Properties of LP 349-25AB 



Property 


Median 


68.3% c.l. 


95.4% C.L 


Tucson 


models CBurrows et al. 1997^ 


LP 349-25AB system 


Age (Gyr) 


0.141 


-0.019,0.023 


-0.04,0.05 


q (Mb/Ma) 


0.863 


-0.019,0.013 


-0.04,0.03 


AT.ff (K) 


137 


-15,21 


-40, 50 


LP 349-25A 


Ma (Mq) 


0.065 


-0.004,0.004 


-0.008,0.008 


Teff,A (K) 


2780 


-30,30 


-50,50 


log(9A) (cgs) 


5.02 


-0.04,0.04 


-0.08,0.08 


Ra (Rq) 


0.130 


-0.003,0.003 


-0.005,0.006 


LiA/Lio'' 


0.75 


-0.52,0.15 


-0.74,0.22 


LP 349-25B 


Mb (Mq) 


0.056 


-0.003,0.004 


-0.006,0.007 


Teff.B (K) 


2640 


-30,30 


-60,60 


log(9B) (cgs) 


5.00 


-0.04,0.04 


-0.08,0.08 


Rb (Re) 


0.124 


-0.003,0.003 


-0.005,0.005 


Lie/Lio'' 


0.96 


-0.08,0.03 


-0.28,0.04 


Lyon models (Dusty; 


Chabrier et al, 


,2000) 


LP 349-25AB system 


Age (Gyr)- 


0.127 


-0.017,0.021 


-0.03,0.04 


q (Mb/Ma) 


0.872 


-0.018,0.014 


-0.04,0.03 


AT.ff (K) 


134 


-15,20 


-40, 50 


LP 349-25A 


Ma (Mq) 


0.064 


-0.004,0.005 


-0.007,0.009 


Teff.A (K) 


2660 


-30,30 


-50,50 


log(9A) (cgs) 


4.93 


-0.04,0.04 


-0.08,0.09 


Ra (Rq) 


0.144 


-0.003,0.003 


-0.006, 0.006 


LiA/Lio" 


0.68 


-0.42,0.24 


-0.64,0.30 


(J - K)a (mag) 


0.795 


-0.002,0.002 


-0.003,0.004 


(H - K)a (mag) 


0.375 


-0.006,0.005 


-0.011,0.011 


(J - H)a (mag) 


0.420 


-0.005,0.004 


-0.008,0.008 


{K - L')a (mag) 


0.544 


-0.009,0.008 


-0.017,0.016 


LP 349-25B 


Mb (Mq) 


0.056 


-0.004,0.004 


-0.007,0.007 


Teff,B (K) 


2520 


-30,30 


-60,60 


log(9B) (cgs) 


4.91 


-0.04,0.04 


-0.08,0.08 


i?B (Re) 


0.137 


-0.003,0.003 


-0.005,0.006 


LiB/Lio'' 


0.95 


-0.18,0.04 


-0.54,0.05 


(J - A")b (mag) 


0.810 


-0.005,0.006 


-0.008,0.013 


{H - K)b (mag) 


0.396 


-0.009,0.006 


-0.013,0.013 


(J - H)b (mag) 


0.415 


-0.004,0.005 


-0.007,0.007 


(A- - L')b (mag) 


0.594 


-0.011,0.012 


-0.021,0.025 



-Abundance of lithium relative to the initial amount (Lio). 
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Table 8. Evolutionary Model-derived Properties of LHS 1901AB 



Property 


Median 


68.3% c.l. 


95.4% c.l. 


Tucson 


models (Bi 


arrows et al. 1997^ 


LHS 1901AB system 


Age (Gyr)- 


0.37 


-0.15,9.63 


-0.21,9.63 


q (Mb/Mj^) 


0.966 


-0.016,0.011 


-0.04,0.02 


AT.ff (K) 


35 


-8,10 


-24, 32 


LHS 1901A 


Ma (Mq) 


0.098 


-0.010,0.004 


-0.018,0.008 


Teff,A (K) 


2960 


-30,30 


-70,60 


log(9A) (cgs) 


5.22 


-0.07,0.03 


-0.14,0.04 


Ra (Re) 


0.128 


-0.004, 0.004 


-0.007,0.009 


LHS 1901B 


Mb (Mq) 


0.095 


-0.011,0.005 


-0.019,0.007 


Teff.B (K) 


2930 


-40,30 


-80,60 


log(9B) (cgs) 


5.23 


-0.08,0.03 


-0.15,0.04 


Rb (Rq) 


0.124 


-0.003,0.004 


-0.006,0.010 


Lyon models (Dusty; 




.2000) 


LHS 1901AB system 


Age (Gyr)^ 


0.28 


-0.08,9.72 


-0.13,9.72 


q (Mb/Ma) 


0.958 


-0.014,0.015 


-0.04,0.03 


AT.ff (K) 


39 


-9,10 


-27,35 


LHS 1901A 


A/a {Mq) 


0.099 


-0.011,0.008 


-0.020,0.013 


Teff,A (K) 


2860 


-50,50 


-90,80 


log(9A) (cgs) 


5.15 


-0.07,0.07 


-0.15,0.08 


Ra {Re) 


0.138 


-0.005,0.006 


-0.008,0.012 


(J - K)a (mag) 


0.787 


-0.001,0.001 


-0.002,0.002 


{H - K)a (mag) 


0.336 


-0.009,0.009 


-0.013,0.018 


(J - H)a (mag) 


0.451 


-0.009,0.008 


-0.017,0.010 


LHS 1901B 


Mb {Mq) 


0.094 


-0.010,0.010 


-0.019,0.014 


Teff,B (K) 


2820 


-40,50 


-90,90 


log(9B) (cgs) 


5.15 


-0.07,0.08 


-0.15,0.09 


Rb {Re) 


0.135 


-0.005,0.006 


-0.008,0.012 


(J - K)b (mag) 


0.788 


-0.001,0.001 


-0.002,0.002 


{H - K)b (mag) 


0.341 


-0.010,0.009 


-0.014,0.019 


(J - H)b (mag) 


0.447 


-0.008,0.009 


-0.016,0.012 



''Both sets of evolutionary models are only computed up to 
an age of 10 Gyr; therefore, this defines the upper limit on the 
model-derived ages. 
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Table 9. Evolutionary Model-derived Properties of Gl 569Bab 



Property 


Median 


68.3% c.l. 


95.4% c.l. 


Tucson models (Bi 


arrows et al. 1997^ 


Gl 569Bab system 


Age (Gyr) 


0.51 


-0.08,0.13 


-0.14,0.38 


q (MBb/MB,) 


0.886 


-0.017,0.021 


-0.03,0.05 


ATeff (K) 


229 


-19,19 


-37,39 


Gl 569Ba 


MBa (M©) 


0.074 


-0.004,0.004 


-0.007,0.007 


Teff.Ba (K) 


2530 


-30,30 


-50,50 


log(gBa) (cgs) 


5.30 


-0.04,0.03 


-0.07,0.06 


flBa (i?©) 


0.101 


-0.002,0.002 


-0.003,0.003 


LlBa/Lio'' 


0.00 


-0.00,0.00 


-0.00,0.00 


Gl 569Bb 


MBb (Mq) 


0.066 


-0.004,0.005 


-0.008,0.010 


T,B,Bb (K) 


2300 


-30,30 


-60,60 


log(5Bb) (cgs) 


5.28 


-0.05,0.05 


-0.09,0.10 


iJBb (Re) 


0.097 


-0.002,0.002 


-0.004,0.004 


LiBb/Lio" 


0.02 


-0.02,0.49 


-0.02,0.82 


Lyon models (Dusty; 




20001 


Gl 569Bab system 


Age (Gyr) 


0.46 


-0.07,0.11 


-0.13,0.33 


q (MBb/MBa) 


0.866 


-0.014,0.019 


-0.03,0.05 


AT.ff (K) 


221 


-18,19 


-36,38 


Gl 569Ba 


MBa (Mq) 


0.075 


-0.004,0.004 


-0.008,0.008 


TefF,Ba (K) 


2430 


-30,30 


-50,50 


log(ffBa) (cgs) 


5.23 


-0.04,0.04 


-0.08,0.06 


^Ba (Re) 


0.110 


-0.002,0.002 


-0.003, 0.005 


LiBa/Lio'" 


0.00 


-0.00,0.00 


-0.00,0.00 


(J - K)b, (mag) 


0.827 


-0.007,0.010 


-0.013,0.021 


(H - K)b^ (mag) 


0.371 


-0.005,0.007 


-0.009,0.018 


(J - H)Ba (mag) 


0.456 


-0.004,0.004 


-0.008,0.009 


(E- - L')Ba (mag) 


0.660 


-0.014,0.014 


-0.03,0.03 


Gl 569Bb 


MBb (Mq) 


0.065 


-0.004,0.005 


-0.008,0.010 


Tcff,Bb (K) 


2210 


-30,30 


-60,60 


log{SBb) (cgs) 


5.20 


-0.05,0.05 


-0.09,0.11 


Rfib (Re) 


0.106 


-0.003,0.002 


-0.005,0.004 


LiBb/Lio'' 


0.00 


-0.00,0.02 


-0.00,0.43 


(J - -ft')Bb (mag) 


0.96 


-0.03,0.03 


-0.05,0.06 


(H - K)Bb (mag) 


0.432 


-0.017,0.014 


-0.03,0.04 


(J - H)Bb (mag) 


0.525 


-0.011,0.013 


-0.02,0.03 


[K - L')Bb (mag) 


0.786 


-0.016,0.015 


-0.03,0.03 



^Abundance of lithium relative to the initial amount (Lio). 



